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ABSTRACT 


Work covered in this report represents a portion of the third year's 
effort of a continuing program tc determine the reliability attributes of 
terrestrial solar cells. Major effort during this reporting period was 
devoted to two tasks: 1) improvement of the electrical measurement instru- 

mentation through the design and construction of a microcomputer controlled 
short interval tester, and 2) better understanding of second quadrant behavior 
by developing a mathematical model relating cell temperature to electrical 
characteristics. In addition, some preliminary work is reported on an inves- 
tigation into color changes observed after stressing. While the accelerated 
stressing of various cell types is continuing, no new results of this activity 
are presented in this report and the reader is referred to earlier reports for 
the latest available documentation- 



SUMMARY 


The third year of the accelerated reliability testing program concen- 
trated on electrical measurement instrumentation and in modeling cell behavior 
in the second quadrant* In addition, some preliminary work was done on corre- 
lating cell color changes with electrical degredatlon. Not reported are 
results of continuing accelerated stress tests on state of the art cells* A 
number of new cells were added to the program, but not in time for sufficient 
data to be obtained, while the older cells are undergoing extended test peri- 
ods and new data is not yet available on them* 

The all-digital, microprocessor controlled, short interval tester, which 
was designed and fabricated on the program, has replaced the manual measure- 
ment procedure formerly used* This has improved measurement accuracy and 
repeatability, reduced measurement time, and through coordinated data manage- 
ment procedures, eliminated data errors. A complete description of the tester 
including schematics and software is given and its operating procedures 
described. 

A computer model, based on the thermal and electrical properties of the 
cells and encapsulating materials, was developed to relate cell temperature to 
electrical characterstics in the second quadrant. This model adequately pre- 
dicted the behavior of both encapsulated and unencapsulated cells, although 
accurate temperature measurements on encapsulated cells were difficult to 
obtain. In addition, only cells of one type were used for comparison and 
other cell types may require different parameter values for fitting. Use of 
the model should permit the prediction of a cell’s sensitivity to degredatlon 
in the second quadrant. The computer program is listed together with a 
description of its operation. 
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In an effort to determine the possibility of quantifying AR coating 
effectiveness by color measurement some preliminary work was done using an IBM 
7400 spectrophotometer system. 
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1.0 INTRODUCTION 


This is the Third Annual Report on the Investigation of Relabillty 
Attributes and Accelerated Stress Factors on Terrestrial Solar Cells, a 
program being conducted by Clemson University for the Low Cost Solar Array 
(LSA) Project of the Jet Propulsion Laboratory. The First Annual Report was 
an introductory report which discussed the philosophy of accelerated testing 
as applied to solar cells, reported the results of applying certain test pro- 
cedures to four types of state of the art solar cells, aud concluded with a 
preliminary reliability qualification test schedule. Also included in the 
first year's program were the development of techniques for reproduclbly 
measuring electrical cell parameters and the ioq>lementation of a computerized 
system for reliability data management and analysis. 

The second year program refined and extended many of the techniques 
developed under the first year's program. Additional accelerated environ- 
mental testing was performed, both on the original four cell types and on 
three additional cell types. These additional cell types Involved advanced 
crystal growing and assembly methods, but did not significantly differ in 
their metalizatlon systems. As a result of this additional testing, several 
changes were made to the preliminary qualification test schedule resulting in 
a reduction in both the number of tests and the number of cells required per 
test. During the second year's effort it became evident that, in addition to 
studying first quadrant failure modes, the contract should spend some time 
characterizing the second quadrant. Results from the field were pointing to 
inadvertent operation in the second quadrant as the major cause of field 
failures. Consequently a second quadrant subtask was initiated. 
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It was also apparent during the second year of the program that It would 
be necessary to test numerous different cell types as the deslgn-co-cost 
philosophy of the LSA program led to different and less expensive metallzatlon 
systems. In order to cope with additional quantities of cells, on top of 
those already In the program, revision of the electrical measurement methods 
was required and a measurement method subtask was Initiated. 

Because of these two factors~second quadrant Induced field failure and 
the need for Improved measurement isethods — most of the effort In the third 
year of the contract, and consequently the bulk of this annual report, has 
been devoted to these two subtasks. The statistical analysis of accelerated 
stress test results, which domlnnted the first two reports, Is thus oiissing 
from this report. However, a comprehensive review of the accelerated testing 
program was presented in a workshop held May I*-2, 1980 at Clemson and 
proceedings are available as a separate JPL document*. Furthermore, a number 
of new cell types with radically different metallzatlon systems were added to 
the program too late for any significant amount of data to be available for 
this report. Table 1.1 Is a summary of the physical characteristics of all 
cell types presently In the program. 

Two additional Items dlsc**^sed In this report are a technique for measur- 
ing the loss of antireflecting coating during stress testing, and an update on 
the data collection system, which has changed somewhat as a result of the 
Improved measurement Instrumentation. 


*JPL Low-Cost Solar Array Project Document 5101-163. 
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2.0 DEVELOPMENT OF A MATHEMATICAL MODEL TO ASSESS THE RELIABILITY 

CONSEQUENCES OF OPERATION IN THE SECOND QUADRANT -- Subtask 1 - 

R. A. Hartman 

2.1 Introduction 

It has been shown, through an extensive laboratory accelerated test pro- 
gram, that the application of high temperature, high humidity, and high rates 
of temperature change, either singly or In combination, can result in irre- 
versible change In the electrical characteristics of solar cells. These 
changes often manifest themselves In a gradual degradation with time of a 
cell's power output. In addition, at sufficient stress levels or if stressed 
for a sufficient length of time, cells may display a catastrophic failure mode 
in which the cell becomes either open circuited or shor circuited. Although 
extrapolation of the observed effects In the laboratory to use conditions Is 
difficult because of uncertainties In the acceleration factors, the best 
estimates of the amount of power reduction which could be expected from state 
of the art cells over a 20-year period gave no cause for alarm. In a number 
of field Instllations , however, a greater than anticlpatd power reduction has 
been observed. Subsequent Investigation of the field performance of deployed 
arrays led to the conclusion that the problem stemmed from cells becoming 
reverse biased and operating in the second, rather than the first, quadrant of 
the VI characteristic. When this occurs the affected cell dissipated an 
Increased amount of power and as a result increases in temperature. An ex- 
amination of operating arrays that were not designed to avoid second quadrant 
operation, using infrared scanning instrumentation, shows dramatically that 
numerous cell or parts of cells are appreciably hotter than their surround- 
ings and consequently are operating In the second quadrant. It can be 
recognized that when this occurs the affected cells are in effect experiencing 
accelerated stressing under use conditions and should experience the same 
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acceleratLu degradation and catastrophic failure rates experienced In the 
laboratory. When a cell open circuits as a result of the high stress level, 
power is of course lost from the entire series string and, depending on the 
particular module and array design, rather rapid and irreversible system power 
loss can result. This section of the report examines factors governing the 
temperature of cells operating in the second quadrant and develops a computer 
model to relate temperature to electrical characteristics. 

A cell will operate in the second quadrant whenever the system passes 
current through the cell in excess of the cell's short circuit current, 

Ig(,. There are a variety of reasons why this might h<«ppen; external 
events, such as bird droppings or fallen leaves may cause partial shading of a 
cell; something could be wrong with the module Itself, such as Interconnect 
failure or a cracked cell; or parametric degradation could occur resulting in 
a change in a cell's electrical characteristic. 

When the operating point shifts from the first to the second quadrant, 
appreciable power is dissipated in the cell. Figure 2.1 Illustrates this 
clearly. It shows the V-I curve of a particular type of solar cell in the 
first and second quadrant. Constant power hyperbola's are superimposed to 
permit estimation of the power which is generated, and therefore must be 
dissipated, at any point along the curve. While this particular solar cell 
can supply about I Watt to an external load when operating in the first 
quadrant, the power that is dissipated in the second quadrant can easily be 10 
to 50 times that amount and in the vertical region of the curve (breakdown) 
will be only limited by the external circuit resistance. Solar cells are 
current sources and the physical configuration of series parallelling will 
determine the maximum current that can occur. (2.1) 
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CELL CURRENT lAl 


The thermal properties of the module determine the shape of the VI char- 
acteristic that will be obtained in the second quadrant. The cell in Figure 
2.1 was suspended in free air without heat sinking. Figure 2.2 shows two 
characteristic traces of a cell. The top trace was obtained as in Figure 2.1 
while the bottom trace was obtained using a pulsed high current fixture and 
curve tracer. The high current fixture permits the solar cell characteristics 
to be taken under pulsed conditions of up to 200 amperes peak current. The 
pulse width is 300 sec at a 60 Hz repetition rate, resulting in a duty cycle 
of less than 2%, and negligible heating. The bottom trace therefore corre- 
sponds to a cell mounted in a module having an infinite heat sink, i.e. a heat 
sink that would not allow a rise in temperature above ambient no matter what 
power was dissipated in the cell. The trace of an actual cell in a module 
will always be above the bottom trace, since the module is not an infinite 
heat sink, and below the top trace because the thermal conductivity of the 
encapsulant is better than the thermal conductivity of air. 

Second quadrant effects were first described by F. A. Blake and K. L. 
Hanson in 1"69. (2.2) The phenomena was "discovered" during studies of prob- 
lems associated with high voltage solar arrays. They concluded that the prin- 
cipal danger of operation in the second quadrant was that it could amplify a 
small and sometimes temporary condition into a complete open circuit failure. 
They further noted that a small number of cells in parallel accentuate the 
severity of second quadrant operation and pointed out the use of by-pass 
diodes. R. N. Diamond and E. D. Steele (2.3) described the possibility of 
integr:’tng a diode for each cell into the back of its neighbor. The article 
discussed the construction in detail and showed V-I curves with and without 
diodes. In the New Technology Section of this report is a description of a 
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CELL VOLTAGE (V) 


Figure 2.2 Comparisons oE Pulsed and Steady State V-I Characteristics 
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further development of the Integrated diode cell. R. N. Diamond (2.4) further 
described radiation effects and diode characteristics. P. L. Jett and J. L. 
Miller (2.5) analyzed three proposed configurations for the orbital workshop 
solar array for hot spot sensitivity. The configurations studied are strings 
of 154 by 1, 2, or 4 cells. Excessive heating was formed when three out of 
four parallel cells were shaded. Further analyses for skylab were made by N. 
B. North and D. F. Baker (2.6). N. R. Garner (2.7) applied a statistical 
analysis technique to determine the maximum number of failure that could be 
expected during a space ship maneuver that casts a shadow on the array. 

H. S. Rauschenbach and E. E. Malden (2.8) discussed breakdown phenomena 
in reverse biased solar cells. They found short circuit failures and occa- 
sional silicon melts. M. Sayed and L. Partaln (2.9) discussed the effect of 
the second quadrant mode of operation on cadlmum sulfide cells. 

R. G. Ross (2.10) described the implications of reverse bias on module 
and array reliability. He pointed out that a proper design of series/parallel 
configurations could reduce the problem. 

Three papers in the German literature by K. Schneider and W. Schultze 
(2.11), by J. Rath and B. Schultz (2.12) and by B. Gohrbrandt and B. Georgens 
(2.13) mainly reviewed the hot spot problem and Its possible solutions. 

The present work, described in this report, develops a mathematical 
model, based on heat flow considerations, which relates cell temperature and 
the shape of the VI characteristic curve In the second quadrant. Such a model 
Is unique and has not been found in the literature. 

2.2 Breakdown Mechanisms 

Figure 2.2 shows a great difference In breakdown voltage when measured 
under pulsed and steady state conditions. The breakdown mechanism associated 
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. pulsed measurement (lower tract), 1- ralanch multiplication. The volt- 
age 'icrops the p-n junction increases ♦ :he point where the electric field in 
the space charge region is sufficie :, t :i create additional carriers by impact 
ionlratioi. This type of junct^ c i :‘iikdow\i has a positive temperature coef- 
ficient, i-e., the breakdown vo . ';;-3 increases as the temperature Increases, 
and is the type of breakdown n tmally found in semiconductor diodes. Two 
other breakdown mechanisms are Zener and thermal. Zener breakdown is caused 
by carriers tunnelling through the barrier and is only found in semiconductor 
devices that have extremely high doping concentrations on both sides of the 
junction. These high levels of doping are not found in silicon solar cells. 
The last of the three breakdown mechanisms involves thermal instability and in 
discrete devices is referred to as second breakdown. Second breakdown is 
often observed in transistors and SCR's, but not in diodes where the pulsed 
and steady state breakdown values are the same because of low leakage currents 
and good heat sinking. 

This is not the case fo*. solar cells, however, where steady state break- 
down is appreciably lower than pulsed breakdown because of the tremendous 
amount of power which can be dissipated in the second quadrant (reverse bias 
p)Olarity) by virtue of the cell being Illuminated. The power dissipated in 
the illuminated cells is roughly the short ciruic current, times the 

voltage across the cell. Since the cell is much hotter in steady state, due 
to the higher dissipation, and yet its breakdown voltage less, the breakdown 
mechanism must not be avalanch multiplication which exhibits a positive tem- 
perature coefficient. 

The ruechanlsm which causes the current to increase as the reverse voltage 
of Figure 2.2 approaches 15 volts is the increase in leakage current due to 
the thermal generation of carriers* 
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At this point the thermally generated leakage current has become an 
appreciable fraction of the total current. As previously discussed, the heat 
which thermally generated the hole-electron pairs comes from t s photon 
current and the reverse potential. In silicon, the thermal generation of 
hole-electron pairs results in an intrinsic concentration of electrons and 
holes per cubic centimeter given by 


n^ = 1.69 X 10 


19 


[3/2 


-0.55/kT 


i300 


10 -3 

At SOC^K the intrinsic concentration is 1.38 x 10 cm , and is negligible 
compared to ordinary doping levels. Consequently any thermally generated 
leakage at this temperature will be orders of magnitude less than the reverse 
diode diffusion current derived from doping levels, which, in turn, is orders 
of magnitude less than the photon generated current. However, n^ increases 
rapidly with temperature, approximately doubling for every 11®C rise in tem- 
perature, and at high temperatures thermal generation will become the dominant 
process of hole-electron production. 

Figure 2.3 is a plot of the temperature at which the concentration due to 
thermal effect, n^^, becomes equal to the original concentration of the base 
as a function of the base doping level. This temperature is termed the 
intrinsic temperature, 0£. It can be seen that for normal base doping 
levels the intrinsic temperature is the order of 300*^0. Below 6^, the 
carrier concentration in silicon is relatively temperature independent, while 
above 0 jj^ the power dissipated in the slice due to the photon and thermally 
generated reverse current causes the cell to heat unifo**mly and operation is 
Inherentlv stable. As the reverse voltage is further increased, the localized 
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temperature at some point will come close to, or exceed 0£- When this 
occurs, thermally generated current becomes appreciable compared to the photon 
current, increasing the localized power dissipation and temperature still 
further, resulting in a positive feedback situation capable of producing 
thermal runaway with an increase in current density by 2 or 3 orders of 
magnitude. It is said at this point that a localized misoplasma has formed. 
Any inhomogeneity that results in a localized temperature in excess of the 
intrinsic value can lead to mesoplasma formation and thermal runaway. This 
inhomogeniety could come about as a material property through structural or 
doping imperfections. Mesoplasma formation will depend on the ratio of the 
rate of power application to the thermal relaxation time of the cell. Thus a 
mesoplasma could be initiated if the current were suddenly Increased to some 
value larger than Ig^, but at the same time, the cell could also operate 
satisfactorily at this same current level had it been approached gradually. 
Destruction of the cell in Figure 2.4 was caused by the sudden application of 
two times Igc> whereas it had operated previously in a stable fashion at 
this same current level for an extended period. An encapsulated cell is less 
susceptible than a bare cell to mesoplasma formation because of the heat 
sinking properties of the encapsulating materials. As a matter of fact, in 
our work thermal runaway and subsequent mesoplasma formation were never 
observed in encapsulated cells, and even for unencapsulated cells it takes 
rigorous application of power to obtain it. The cell in Figure 2.4 had a 
voltage in excess of 60 volts suddenly applied and it took appoximately 3 
seconds before the mesoplasma formed and melted the hole. Proper design 
should be capable of avoiding the generation of such sudden high voltages and 
currents within a module. 
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Figure 2.4 Hole Formed in Cell as a Result of Mesoplasma Formation 
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On the other hand» long before mesoplasma formation *iith its thermal run** 
away is reached on the VI charactelstlc» non uniform heating and hot spot 
formation can occur. Hot spots will occur as temperatures approach but 
where the carrier concentration does not necessarily have the strong exp>nen- 
tial variation with temperature that it does above £• Hence, a hot spot is 
inherently stable. It is in this sub^lntrinsic region that most of the 
reverse voltage reliability problems occur as a result of localized cell tern** 
pc natures of 100-200*C. The higher the temperature the stronger the carrier 
concentration dependence on temperature and the more localized the heating. 

Because second quadrant breakdown is a thermal effect, it is Important to 
understand the temperature dependence of the reverse current. To this end it 
can be assutned that the reverse current consists of several independent compo** 
nents, each with its own specific temperature depend^.nce. The assumption of 
independence will be valid except near avalanch breakdown where multiplication 
should be taken into account. 

The largest current component » at least at low voltages, is the photon 
generated current which has a positive temperature coefficient due to narrow- 
ing of the bandgap at elevated temperatures. When the bandgap becomes smaller 
the absorpt .on edge shifts to longer wavelengths giving rise to Increased 
current. In the temperature range of Interest the effect will be linear 
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where Is sc room temperature. The actual temperature coefficient 
depends on Che spectrum and Is less than H of I sc per ^C. 

Another temperature dependent component of current la the diode satura- 
tion current, also kno%ni as Che diffusion current. For a one-sided p'^n 
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abrupt junction^ similar to a shallow, heavily diffused junction, the current 
density is given by 




where the quantities D^, p^^, and * ^^p^p^ temperature dependent. 

Assuming, that proportional to some power of the temperature such that 

n /i \ X where s • constant, then, 

P P 



exp 
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and the predominant temperature dependence will be that of the exponential 
term and a plot of the logarithm of current vs. l/T should be a straight line 
with slope ■ -Eg/k. 

The third temperature dependent component of current is the diode recom- 
bination-generation current. The temperature dependence of this component can 
be discussed in terms of a single trap located in the middle of the energy 
gap. For this case it can be shown that the rate of generation of electron 
hole pairs is 
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the carrier thermal velocity, is the trap density, the trap energy 
level, the intrinsic Fermi level, and n, the intrinsic carrier density. 
Assuming that "v:j and that the trap energy level is located at the 

intrinsic Fermi level then, 
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The recombination-generation current density is proportional to the space 
charge depletion width, W, so that 

•^rg ^ 

Thus a plot of the logarithm of this component vs. 1/T will also be a straight 
line^ but with a slope » -Eg/2k. Since W is a function of the square root of 

voltage it should be noted that for an abrupt junction, 

J -V V" 

rg 

2.3 Basic Equations 

In order to derive an equation that describes the heat f)ow in a solar 
cell in the second quadrant it is necessary to make some a priori assumptions. 
First it is assumed that the cell Is r^und and that the temperature distribu- 
tion has radial symmetry. This means mathematically that if a hot spot 
exists, it will be located in the center of the cell. Secondly, it is assumed 
Chet there is no temperature difference between the top and the bottom of the 
cell, i.e. all vertical temperature gradients are zero. These two assumptions 
result in a one dimensional time dependent problem. Figure 2.5 indicates the 
assumed cell geometry. Equations may now be set up governing heat flow into 
and out of the differencial ring element, as well as the heat generated within 
Che ring as a result of current flow; 

Heat conduction into the ring from the center, 

q(in at r) = -2^rwic — . 

Heat conducted out of Che ring toward Che outside, 

q(out at r + Ar) « 2i (r + Ar)w < —\ 

^^'r + A r 

Heat radiated to the ring from Che light source, 
ql^radiaced) ■ 2^rULL)Ar 
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Heat lost by the ring by conduction to ambient air, 


q(lost) » -2TTrh . 

Heat generated by current flowing within the ring, 

q(current) * 2:TrjVAr 

Heat stored in the ring, 

3T 

q(stored) * 27rrArwpc — 

vJ L 

where w is the thickness of cell (cm), 

< is the thermal conductivity of the silicon (W/cm®C), 
j is the current density of the current across the Junction (A/cm^)^ 

V is the applied voltage (V), and 

h is the heat exchange between cell surface and surroundings (W/cm2®C) 

ILL is the incident energy flux from the sun (W/cm^®C) 
p is the specific mass of silicon (g/cm^) 
c is the specific heat of the silicon (Ws/®Cg). 

Combining all these terms into a single eq\iation ,<»lves, 

/t ^ 1 3T ^ I 3k ^ jV - h(T - T . ) + ILL I 3T 

\ / sW 

where ^ ^ » the thermal diffusivity. 

►'C 

This equaclo.i decribes the behavior of the unencapsulated solar cell at 
0 < r < R, 

To include the effects of a layer of encapsulant the mode) Is changed 
slightly by exchanging heat between the cell and the encapsulanr Instead ^f 
between the cell and the surroundings as shown In Figure 2.6. An additional 
set of equations is added that describes the heat transport In the encapsulant 
and between the encapsulant and the surroundings- These changes result In the 
following set of equations: 




For the silicon slice: 
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q(out at r + Ar) = -2:r(r + Ar)w k ^ 
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For the encapsulant: 
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q(out at r + ^r) = -2Tr(r^ Ar) w < 
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where the subscript s denotes the silicon cell and e the encapsulant » 
and hgg stand for the heat exchange between the cell and the encapsulant 
and between the encapsulant and the air respectively. Setting up the heat 
balance as before gives: 
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Note that the square of the first derivative does not occur In the encapsulant 
equation, since the thermal conductivity was assumed Independent of tempera- 
ture. For most encapsulant materials this seems to be justifiable. In 
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addition, the illumination term has been moved to the encapsulant equation to 
permit calculation of temperatures where there is encapsulant only, such as in 
the area directly around the cell that is studied. This might not be justi- 
fied in high density modules. 

The model does not allow for a variation of applied voltage across the 
cell due to surface resistivity. Actual variations of the order of 5% were 
measured between the voltage at the voltage probes, and that directly across a 
hot spot. This effect was neglected, since introduction of it would compli- 
cate the computations substantially. At high currents the resistive drop can 
shift a hot spot towards the current probe at the back contact. 

The next section discusses some of the variables in the equation, like 
the current density and thermal conductivity, that are strong functions of 
temperature. 

2.4 Temperature Dependent “Constants” 

Several important ‘’constants” in the equations derived in the last 
chapter such as the current density and the thermal conductivity, are really 
strongly temperature dependent. 

Figure 2.7 shows the temperature dependence of the I-V plot. Clearly the 
current in the second quadrant is a strong function of temperature. These 
curves were obtained with the curve tracer and pulsed high current fixture 
while the cell was mounted on a hot plate equiped with vacuum hold down, 
voltage and current probes, and a thermocouple touching the back of the cell. 
The illumination was lOOmW/cm^ from four ELH lamps. Figure 2.8 shows the 
construction details of the fixture used. 

If in figure 2.7 a constant voltage is selected then it is possible to 
obtain plots of current versus temperature as si ^wn in Figure 2.9 for a 1 volt 
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Characteristics as a Function ot Temperature 
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reverse bias. Appropriate curve fitting techniques applied to a particular 
cell's data plotted in the form of Figure 2.7 and 2.9 led to an emperical 
equation for the current density of: 

J = 1.933 X lO'^ + 8.127 x 10~^T + 1.453 x 10“^exp(4. 585 x 10"^T) 

+ 2.297 X 10"\ exp(1.208 x 10~^T) 

This expression indicates that for a given value of voltage there are three 
regions having different temperature dependence — a linear region and two 
regions of different exponential dependence. These three regions can be 
clearly distinguished in Figure 2.9. 

Region A (T < 180®C): In this region the linear component of the current 

density dominates. The coefficient shown is O.A%/®C which is of the order of 
the previously estimated 1%/®C. 

Region B (180®C < T < 250®C): In this region the reverse current gener- 

ated by traps is dominant and, as shown earlier, should have an exponential 
temperature dependence of -Eg/2k. As shown in Figure 2.9 the experimental 
data in this region displays this type of behavior exactly. It was also shown 
theoretically that the recombination-generation current for an abrupt 
junction, such as farmed by a very thin diffused layer of high conductivity, 
should exhibit a voltage dependence. No explanation can be given at this 

time for the apparent linear voltage dependence observed. 

Region C (T > 250®C): In this region saturation dominates. The lower 

boundary coincides with the intrinsic temperature for 10 ohm-cm mrterial (base 
region). In this region the exponential slope approximate “Eg/k as was 
theoretically predicted. 
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Avalanch multiplication effects were not taken into acccunt when fitting 


curves of Figure 2.7, i.e. only the linear portion of these curves was 
used. This will affect the validity of the model near breakdown. However, in 
all the simulation runs to be described, voltages were kept well below 
breakdown. The breakdown voltage for this particular cell type was around 40 
volts with a positive temperature coefficient of around 0.01 V/®C. 

Values for the temperature dependence of the thermal conductivity can be 

obtained from the literature. In particular, E. C. Haesell (2.14) suggested 

the form k* - k /(T + T ) with < * 250 W/cra and T » 205"*C and this was used 

o o o ® 

in the simulations. 

Other factors such as the specific heat, are not really Independent of 
temperature. However these have a much weaker temperature dependence than the 
current density or the thermal conductivity in the region of interest. The 
model was only used to calculate steady state temperature* It might be 
desireable to include temperature dependence if the model were to be used for 
tranisent calculations. It should also be pointed out that the heat exchange 
between a solid surface and its surroundings is most likely not linear 
process and may well be a function of factors other than temperature. The 
assumption of a linear dependent heat loss mechanism is an a priori 
assumption. 

2 . 5 Boundary Conditions 

’"or the unencapsulatod cell boundary, condltlors must be defined at the 
center and the edge of the cell. For the encapsulated cell boundary, condi- 
tions must be defined at these points and at the edge of the encapsulant. 

Because rotational symmetry was assumed, the first derivative, , must 
equal zero at the center of the cell. This causes the deletion of the two 
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terms containing this derivative in the cell equation and of one term in the 


encapsulant equation. The rotational symmetry is a geat computational 
advantage as will be obvious after reviewing the Thomas algorithm, described 
in Appendix A. 

Using the methods described in Section 2.3 Equation 2.2 may be evaluated 

at the boundary between the cell edge and the ambient or encapsulant. After 

some manipulation the following expression is obtained for the case of the 

encapsulated cell: 
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The equation for the unencapsulated cell is analogous. The third therm on the 
left is the amount of heat that leaves the cell v^a the cell edge. This 
amount is added to the next interval, the arrow in Figure 2.10 symbolizes this 
flow of heat. 

The final boundary condition that has to be set is the end point of the 
encapsulant to be simulated. An extra ring is added to the encapsulant to 
indicate a continuation of it. Temperature was first extrapolated in the 
following fashion. 


T = ?T - T 

*N+1 N-1 , 


but this tended to introduce an oscillation at the edge of the modeled 
encapsulant. An obviously less accurate, but more stable solution was 
obtained by using 
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.10 Encapsulated Cell Space Grid 



In order to fit the results of the analysis to the measured temperatures 
and I-V curves, it was necessary to introduce a perturbation in the current 
density. This was done by Introducing a Gaussian spatially dependent term in 
the exponent of the leakage current term of the form: 

8 + Y exp 

The spread, a , of the Gaussian was estimated from experiments to be one third 
of the cell radius. 6 and y are arbitrary constants which can be used to fit 
the model to the results. The term 8, if decreased, results in a cell that is 
less leaky that the typical cell, and consequently a cell with a higher base 
doping, as explained in Section 2.2. The factor y was varied between 0 and 
0.6, with a typical fit being obtained for a value of Y « 0.4. If Y is 
increased, a localized spot that is lighter doped, or for some other reason 
leakier, is Introduced at the center of the cell. 

Because of the assumption of circular symmetry the location of the leaky 
spot (defective area) is restricted to the center of the cell. This, of 
course, is not what is observed experimentally, but it is necessary in order 
to solve the equations. There is little loss of generalization, however, as a 
consequence of this assumption, mainly because the exact location of the hot 
spot does not have much effect on the cell's VI characteristic. The only 
possible exception might be a hot spot which formed at the cell edge. The 
model for encapsulated cells consists of two linear equations, and calculation 
of one I--V curve takes 7.5 minutes of CPU time using an IB 370/3033 computer. 




Locating the hot spot anywhere else on the cell would take ten hours of CPU 
time by conservative estimate! In a similar fashion the model is restricted 
to circular cells, but this is presently the most popular form factor. As 
square cells proliferate their second quadrant performance will need to be 
stud led. 

The computer programs for modelling both unencapsulated cells and 
encapsulated cells in the second quadrant are listed in Appendix 0. 

2*6 Unencapsulated Cell Results 

Measurement Methods — Figure 2.11 shows four experimentally measured I-V 
curves of cells in the second quadrant. These- curves were obtained with an 
xy-recorder, a programmable power supply and a light source. The light source 
was a single ELH lamp overhead projector which gave an Illumination level at 
the cell of 100 mW/cm2, The power supply was programmed as a current supply 
and could be varied, either through the use of a variable resistor, or by 
means of a function generator. The slowest voltage ramp of the function 
generator was 5mV/sec, which was too fast to generate a steady state V-I 
curve, and consequently this method was used only for demonstrations. 

Figure 2.12 schematically illustrates the method used. The front of the 
cell was placed junction side down on the 100mW/cm2 light source. The back 
of the cell was covered with a temperature sensitive phosphor which gave an 
indication of the temperature distribution. The phosphor used was #2090, 
fabricated by the U.S. Radium Corporation. Appendix C of the Second Annual 
Report contains information on temperature sensitive phosphors. These 
phosphors are not readily available in small quantities, but sample quantities 
for evaluation purposes may be obtained from the U.S. Radium Corporation. 

Figures 2.13. (a), (b), (c), and (d) illustrates the use of phosphor 
coatings to detect the development of a hot spot. Figure (a) shows the cell 
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Figure 2.12 Schematic of Second Quadrant Measurement Method 
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at zero voltge, Figure (b) just before the knee in the second quadrant, 

Figure (c) just above the knee at moderate currents, and Figure (d) at high 
currents. 

In some cases the peak temperature on the back of the cell was estimated 
using a wax temperature indicator, such as Omega's melting point standards. 
Good correlation was obtained between the two methods. Temperature measure- 
ments with melting point standards are very accurate, but are only available 
in intervals. Normally during measurement, a thermocouple was attached 

to the rim of the cell. Direct measurement of the hottest spot using a 
thermocouple was attempted, but was not successful. The problem was in 
maintaining good thermal contact for a prolonged period without cooling the 
hot spot by conduction along th® thermocouple. It is possible to measure the 
temperature remotely using an infrared scanner as described by C. A. Lidback 
(2.15). There are, however, some problems with this technique. Tne lead on 
the back of the cell has emmissivity of only approximately 10%, so it is 
necessary to paint the back of the cell with a paint having better emmissive 
properties. Figure 2.14 shows a thermograph obtained in this manner.* The 
method is reasonably accurate (*flO'^C) and at the same time supplies a 
pictorial description of temperature variations across the cell in terms of 
graph tones. Drawbacks are that the positioning of this particular scanner 
over the light table with the cell on it is cumbersome, since both the light 
table and the scanner are bulky objects. 

Model Prediction — As has been pointed out the model permits calculation 
of the IV characteristic curve. For unencapsulated cells the IV curve was 


* The IR scanner used for this work was a model 900 manufactured by UTI and 
owned by North Carolina State University, Raleigh, NC. 
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Figure 2.14 IR Scanner Thermograph of an Unencapsulaced Cell 
in the Second Quadrant 
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Figure 2.15 Calculated VI Characteristics of an 
Unencapsulated Cell for y* 0.2 
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2.16 Calculated VI Characteristics of an 
Unencapsulated Cell for '!■ 0.4 
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obtained from 40 points calculated at O.IA intervals starting at a current of 
1.3A. Figure 2.15 shows IV curves for which the amplitude of the Gaussian 
perturbation, Y , is 0.2, while the temperature dependent parameter of the 
leakage current, 3 , is varied between 0.7 and 1.2 in intcrrvals of 0.1. Figure 
2.16 shows the same but for Y * 0.4. Trace 1 in Figure 2.11 is from the cell 
used to model the current as function of temperature and voltage (Cell 1), so 
this curve should closely ma^'ch a c. ’“ve with 6 * 1.0. Figure 2.17 shows the 
measured curve of Cell 1 superimposed on the calculated curves for 3*1. For 
Y » 0.4 there is satisfactory fit between the calculated curve and the 
measured curve up to .5A, aporoximately two times Figure 2.18 shows 

a dia;^ram of Cell 1 with a description of the visual observations made as the 
current through the C-11 was increased. Under the UV lighting conditions 
used, the transition point for the phosphor was approximately 200°C. In 
Figure 2.18 this transition is indicated as the appearance of a hot spot. The 
coeyistance of two hot spots at moderate current levels, followed by the 
disappearance of one at higher currents, is often observed. A shift cowards 
the back currenc probe is not uncommon at high current loads, and probably is 
due to a resistive voltage drop in the base ind back metalizations of the 
solar cell. Measurements of the resistive voltage drop gives va^ue in the 
orcar of .2 - .3 V, an effect that was not incorporated in the model, as 
mentioned earlier. 

Table 2.1 compares measured and calculated temperatures. The agreement 
between the model and the measured values is remarkable, considering this was 
not a cel] with a hot spot centered in the middle of the cell, as was assumed 
in the model. 

Figure 2.19 shows the calculated temperature pro les along a radius, for 
::)» 1 and 'f * 0.4, at different current levels. The eflect at the edge of the 
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Figure 2.17 VI Curve of Cell 1 Compared wich 
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Table 2.2 Measured and Calculated Temperature at 
the Edge of an Unencapsulated Cell 
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cell that had been observed experimentally can be seen in Table 2,2. The 
higher temperature observed at lower currents can be attributed to secondary 
hot spot formation in this region. The lower temperatures at higher currents, 
on the other hand, are probably caused by the geometry of the hot spot and how 
it is affected by Che thermocouple. 

Comparing the VI curves of Figure 2.15 and 2.16 with the remaining three 
traces of Figure 2.11 permits an estimation of the parameters for these 
curves. A reasonable fit can be obtained for cell 2 at 3 ■ 1.1, ^ ■ 0.2 and 
for cell 3 and 4 at * 0.8, Y - 0.3. Cell 3 shorted at the end point of the 
trace, and as far as could be determined this was caused by molten lead at the 
table attachment point. For cell 4 the measured temperature at 1.58A was 
between 206®C and 220®C, while the calculated center temperature at 1.5A was 
190°C and at 1.6A was 215°C. At 1.76A the measured temperature is smaller 
than 258®C while Che calculated temperatures at 1.7A and 1.8A are respectively 
234°C and 250°C. Both these observations indicated that the current function 
could be modified Co provide a reasonable fit of the data. 

Parametric Changes — It is possible to vary a number of parameters 
within the model. For example, variations could be made in light intensity, 
partial shading (radial symmetry only), temperature, and cell thickness. 

Figure 2.20 shows the influence of varying cell thickness on the VI cu ve. 
Thicknesses of 0.02, 0.04 and 0.08 ram were chosen. As the cell thickness 
becomes smaller, the lateral transport of heat becomes smaller, causing a 
lower breakdown voltage. The calculated temperature profiles at 2A, shown in 
Figure 2.21, illustate this. Clearly the temperature distribution becomes 
Less uniform at smaller cell thickness. Experiments with a cell, which was 
artificially chickened by attaching a thin brass foil to its back, confirmed 
this. Breakdown voltage in this case was increased from 15V to 21V. 
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Figure 2.20 Calculated VI Curves for Different Cell Thicknesses 
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Figure 2.21 Calculated Temperature Profile for Different Cell Thicknesses 



2.7 Encapsulated Cell Results 


Construction — Figure 2.22 illustrates the construction of the encapsu- 
lated minimodules used in the investigation. These minltnodules were made by 
Springborn Laboratories, Inc., for the Jet Propulsion Laboratory LSA program* 
The cells were manufactured by Solar Power Corp. for the JPL Block I procure- 
ment, and were of the same type that had been studied earlier by Clemson when 
Investigating accelerated reliability testing in the first quadrant. 

The glass superstrata minlraodule consisted of the glass superstrata, a 
layer of clear ethylene vinyl acetate (EVA), the solar cell, a layer of EVA 
containing white pigment with an unwoven glass matt and covered by alumlnimum 
foil. The glass matt is to ensure insulation between the cell and the 
aluminum foil and to add to the mechanical strength. 

The minimodule with the masonite substrate consisted of a film of korad, 
a layer EVA, the cell, EVA with pigment and glass matt, a Super Dorux 
substrate, and an EVA layer. The korad film is a UV screen while the bottom 
EVA layer is to protect the masonite from the elements. 

No data is available on the tempe. ture variation of the EVA thermal 
conductivity. Most probably there is a change in the thermal conductivity at 
higher temperatures when the EVA softens and actually becomes liquid aroond 
105*"C. Another transition point is around 350®C when decomposition starts. 
Figure 2.23 shows a partial trace made with a therraogravimetric analyzer 
(TGA)* of a sample of EVA. Weight loss of the EVA began around 120®C where 
the sample probably started to loose water. At 225®C all the water was gone. 
Decomposition began around 275®C and the sample was completely decomposed at 
550®C. An additional temperature effect is that the EVA becomes less clear 

* Mod^l 951 Thermogra^Timetric analyzer (Dupont) 
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dua to carbonizing, causing a loss of transparency and a drop in the short 
circuit current. These effects were not taken into account in the model, 
however. The model was run primarily to solve steady state problems. In 
these cases the thermal dif fusivities of the silicon and the encapsulant were 
set to fixed values to ensure fast convergence. This makes the value of the 
time increments uncertain, but for steady state problems this is not important 
since the solution is taken at t » 

2.8 Measurements and Results from the Model 

Figure 2.24 and 2.25 show Vl-curves from the glass encapsulated cells and 
the masonite encapsulated cells respectively. In order to relate these curves 
to cell temperature, as was done for the unencapsulated cells, it is necessary 
to have an accurate method of measuring the cell temperature. Unfortunately 
this could not be easily accomplished. In some cases a thermocouple was 
mounted in the center of the masonite substrate touching the back of the cell. 
This procedure, however, infuenced the gas containment capabilities of the 
encapsulant and perturbed the experiments. Some thermocouple temperature 
measurement of the front and back surfaces were attempted, but they only gave 
a qualitative indication of the cell temperature. The development of more 
accurate temperature measuring techniques, such as using an array of very low 
mass thermocouples encapsulated with the cell, was not attempted. An infrared 
scanner is not a particularly useful tool for measuring the temperature of 
encapsulated cells since it measures the surface temperature of the encapsu- 
lant rather than the cell itself. Furthermore, as described in Section 2.6, 
the accuracy of the instrument depends on the surface emissivity and it would 
be necessary to paint the cell with black paint. An infrared scanner, how- 
ever, is a powerful tool for locating reverse biased solar cells in arrays 
deployed in the field. Phosphor decoration also measures only encapsulant 
surface temperatures. The best indicator of cell temperature was thu melting 



of the solder metalization which could be observed at approximately 170®C. 

This effect provided a temperature check which allowed correlation with 
unencapsulated cell results. A total of four miniraociules of each construction 
type were tested in the second quadrant. Two of the glass superstrate type 
(A, B) and three of the masonite substrate type (C, D, E) are shown in 
Figures 2.24 and 2.25. Two additional cells having glass superstrates shorted 
during testing. The first shorted when the llghtsource malfunctioned while 
the cell was biased at 15.7 V and 1.366 A. The generated heat was si^^ficient 
to shatter the top plate of the overhead projector light source. The shock of 
the explosion probably caused the short to occur in the mint module. The 
second cell had a manufacturing defect which caused the juntion to short at 
approximately 7 volts reverse bias. The other masonite substrate mir > .odule 
which was tested, but not shown in Figure 2.25, was accidentally broken. 

The voltage on Cell A (Figure 2.24) was gradually increased from zero. 

At 10.7 volts reverse bias the superstrate cracked due to thermal stress. 

These cracks propagated with time under continued stress. Figure 2.26 shows 
the minimodule two hours after the cracks first started to form. At 15 V 
significant delamination had occurred at the back of the cell and a gas filled 
bubble had formed. At 22 V the cell shorted, probably at a crack location. 
Jelaraination of the top EVA coating can also be seen in Figure 2.26. Lead 
melting was observed at approximtely 18 volts. 

The glass superstrate of Cell B did not crack when subjected to the same 
reverse bias stressing. Instead the cell was pushed downward at a reverse 
bias of 21 volts by gas that formed between the cell and the glass. At this 
point the solder coating was molten. At 18 V a bubble had formed in the back 
EVA. The glass surface did not exceed 149®C during »e test. Figure 2.27(a) 
and 2.27(b) show the onset of gas formation at the top of the cell. The 
photographs were taken a few seconds apart. 


□ GLASS BROKEN 
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Figure :,:4 Measured 71 Curve of Class incapsulated Cells 
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Figure 2.25 Measured VI Curve of Masonite Encapsulated Cells 
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Figure 2.26 Minimodule with Glass Superstate Cracked 
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Masonite substrate Cell C showed degradation of the back EVA at 7V 
reverse bias. A thermocouple » which had been inserted through a small hole in 
the masonite and was touching the back of the cell, registered 120®C. At 
9.4 V the thermocouple Indicated 146**C and delamination had occurred between 
the ^'op of the cell and the top EVA coating as a result of the generation of 
gas. This resulted In the permanent deformation of the korad. Figure 2.28 
shows the appearance of the cell after it had cooled to room temperature. 

Since the bubble clearly was a failure of the encapsulant, the module was not 
subjected to further stesses. 

Masonite substrate Cell D was "llfetested” at 10 V reverse bias for 3 
hours. Figure 2.29(a) shows the severely degraded back coating of EVA after 
this stress. The cell was cracked across the middle as can be seen from the 
front view in Figure 2.29(b). A possible explanation is that water evolving 
from the masonite caused the cell to crack. 

Masonite substrate Cell E was the only encapsulated cell that remained 
reasonably functional after passing the knee of the IV-curve. The IV-curve 
was obtained by stepping in voltage rather than current, but with the power 
supply limited to 3A maximum current. No other ad Jus its were made in the 

power supply settings. At 10 V degradation of the back J7A was again 
observed. At 12 V the cell cracked and melting of the solder was observed 
(170®C). At 16 V small bubbles at the end of the crack appeared. At 18 V 
additional bubbles formed completely around the rim. A small bubble that had 
been formed during encapsulation got bigger, and a strong burned wood odor 
became noticeable. A hot spot developed at just over 21 V. At this point the 
smaller part of the cracked cell probably became Internally disconnected, 
because it appeared noticeable cooler than the larger part. The EVA was 
observed to turn slightly brown all over the cell and dark brown at the hot 

spot. The EVA seemed to flow away from the hot spot and the test was 
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Figure 2-28 Cell C After Second Quadrant Operation 







ter :lnated for fear of fire. Figure 2.30(a) shows the front of the 
minioioduley with the charred hot spot clearly visible and Figure 2.30(b) shows 
the severely degraded back. Some EVA actually flowed from the minimodule to 
the support underneath. 

Since only one cell survived beyond the knee of the IV-curve it was not 
easy to fit the model to the observations. In addition, the end points of 
many traces did not represent true steady state conditions since they were not 
in thermal equilibrium beyond breakdown. Furthermore, many cells were physi* 
cally broken near the end point — a condition that would give a discontinuity 
in the thermal conductivity. 

Figure 2.31 shows calculated IV curves for the glass superstrate module 
for Y » 0.2 for different /alues of 6. Figure 2.32 shows the same for Y • 

O.A. Figure 2.33 and 2.34 show the same curve for the masonite substrate 
module, again for Y - 0.2 and Y » 0.4. 

A reasonable fit for Cell A is 6 » 0.8, Y » 0.2, for Cell B: 3 » 0.7, 

Y ■ 0.2, for Cell C: 1.2, y =0.5 and for E 8 * 0.7, y ■ 0.2. Cell D would 

fit between 3 ■ 0.7 to 3 - 1.0 and Y * 0.2 to y • 0.4. Table 2.3 compares the 

estimated temperatures with the nearest calculated ones. Agreement between 
the calculated and observed values was reasonably good, but not as good as for 

the unencapsulated cells. This can be attributed to a combination of 

measurement difficulties and modeling innaccuracies . 

It is interesting to compare values calculated by the Clemson model with 

data supplied by the Jet Propulsion Laboratory at the Encapsulation Workshop 

held 9/23/80 in Pasadena. Figure 2.35 compares JPL calculations with values 

compiled from the model using glass encapsulation and 3 ■ 1.0, y ■ 0.4. The 

two curves are in excellent agreement over the full range of JPL data. 

Extending the calculations to higher values of P/Pw*^* the Clemson model shows 
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Table 2.3 Estlmas ed Temperatures at Observation Points Compared 
with Calculated Temperatures for Encapsulated Cells 
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Figure 2.31 Calculated VI Char acteri^st ic of a 
Class Encapsulated Cell for **» 0.2 
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figure :.35 Comparison of Clemson '‘oii-l with J?l Prediocion 


a bending back of the curve, indicating that although the hot spot temperature 
is continuing to increase the power into the hot spot is decreasing. This 
thermal runaway condition results as the silicon at the hot spot location 
approaches intrinsic at 25C-300®C. Doubling back of Figure 2.35 occurs as the 
knee of the VI characteristic is approached. 
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3.0 DEVELOPMENT OF IMPROVED SOLAR CELL MEASUREMENT METHODS - Subtask 2 - 
C. R- Saylor, J. F. Christ 

3 . 1 Introduction 

As the number of samples in the accelerated test program increased it 
became quickly evident that the electrical measurement procedure was the key 
to obtaining meaningful results* The measurements made before and after 
stressing had to be highly repeatable in order to detect the small changes 
which occurred as a result of stress* Furthermore, in order to perform the 
quantity of tests necessary to obtain statistically significant results, 
measurement had to be rapid* The potentiometric measurement techniques, which 
were used during the first two years of the program, and which were described 
in the First and Second Annual Contract Reports, were satisfactory as long as 
the volume of cells being tested was small and provided great care was used by 
project personnel during measurement. In order tc obtain high repeatibility 
with a high throughput and inexperienced operators new measurement 
instrumentation was required. Since inexpensive Instrumentation having these 
characteristics was not commercially available a subtask was initiated for its 
design and fabrication and the resulting tester is described in this section. 
Although originally developed for reliability work, the system is flexible and 
by virtue of its low cost, high repeatability, and high throughput should be 
useful in routine industrial measurement applications as well* 

The measurement philosophy of the new instrumentation was to electrically 
characterize the cells rapidly enough that their temperature would not change 
during measurement. In practice this meant making the measureuents in less 
than 1 second* In order to obtain satisfactory light source stability and 
spatial uniformity it was decided to use continuous ELH lamp illumination, but 
with a shutter to avoid cell heating due to illumination. To obtain the 
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desired repeatability and speed an all digital microcomputer controlled system 
was called for. For these reasons the tester may be termed a microcomputer 
controlled, constant temperature, short interval IV tester. The term ’’short 
interval” is meant to differentiate a shuttered system of this type from 
either a flash type tester or a steady state type tester. This approach to 
electrical measurement allows the cell temperature to be controlled prior to 
illumination by air flow rather than the massive water cooled heat sinks 
required previously in the steady state method. The low thermal mass of an 
individual cell should respond quickly to a regulated air stream. Finally, an 
all digital system would eliminate much human interaction and allow the data 
to be introduced directly to the IBM 370 system, where it could be 
statistically manipulated. These general philosophical guidelines led to the 
overall system design of Figure 3.1 

When the shutter is opened a program stored in the microcomputer steps 
the cell through the I-V curve, storing the voltage and current values of each 
data point in memory. Operation is as follows; The desired current value 
through the solar cell is computed in the microcomputer in digital form and 
converted to an analog voltage using the digital controller, which in turn 
controls the programmable power supply. The programmable powei supply elec- 
tronically loads the cell to the current value given. The voltage across the 
cell is measured using an analog-to-digl tal converter whose output the micro- 
computer reads and stores. The process of electronically loading the cell to 
the given current value, measuring the voltage and storing the data, is 
repeated until the entire I-V curve is measured. The T-V curve data points 
that are stored in memory locations can be displayed on the oscilloscope by 
")nnecting two digital-to-analog converters to the X wind Y inputs. Hard copy 
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Figure 3.1 Measurement System Block Diagram 


79 








of the digitized I-V curve can also be plotted using the X-Y plotter. The 
individual data points and/or parameric data can be printed using the teletype 
and also be output on paper tape to be fed into a large digital computer for 
further analysis. Implementation of the system concept of Figure 3.1 into 
specific hardware will now be discussed. 

3.2 System Hardware 

The system hardware can be divided into four major catagories: micro- 

computer, interface hardware, solar simulator and cell holder. Each of these 
will be discussed in detail. 

3.2.1 Microcomputer — The microcomputer used in the short interval 
tester is a single board computer, the SBC 80/10A, manufactured by the Intel 
Corporation (3.1). It is a complete computer system on a single 6.75 x 12 
inch printed circuit card. As can be seen in Figure 3.2, the system is 
designed around three system busses which provide all necess iry communication 
between different parts of the system. 

Intel’s single chip 8-bit n-channel MOS 8080A CPU is the central pro- 
cessor for the computer. The 8080A contains six general purpose registers and 
an accumulator. The six general purpose registers may be addressed individu- 
ally or in pairs providing both single and double precision operators. The 
lb-line address bus allows addressing of up to 64k bytes of memory. An 
external stack, located within any portion of memory, may be used as a last-in 
first-out stack, to provide subroutine nesting capability. A sixteen line 
address bus and eight line bidirectional data bus are used as the interface to 
memory and input/output (I/O) devices. 

The SBC 80/10A contains 48 prograraraab] e parallel I/O lines implemented 
using two Intel 8255 Programmable Peripheral Interface integrated circuit 
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chips. The 3255* s can be programmed in software in combinations of 
unidirectional /bidirectional, input /out put ports. 

The SBC 80/10A contains Ik (1024) 8-bit words of Random Access Memory 
(RAM) and sockets for up to 8k bytes of Read Only Memory (ROM). 

A programmable serial communications interface using Intel's 8251 
Universal Synchronous/Asynchronus Receiver /Transmitter (USART) is contained on 
the board. The USART can be programmed in software to interface to a number 
of communications devices, by programming transmission rates, control 
character format, data format, etc. A jumper selectable 20 ma current loop 
for interface to a TTY is also provided. 

Both the 8251 USART and the two 8255 I/O interface chips are programmable 
and must be configured under software control before they are used. The 8251 
is initialized by sending a mode word followed by a control word to the 
USART's control port. In the short interval tester, the 8251 is configured as 
shown in Figure 3.3 and 3.4. The 8255*s are configured by sending a single 
command word to each of the I/O interface chips* control port and are 
configured as shown in Figure 3.5 and 3.6. 

The measurement system required 17k bytes of RAM, and 9k bytes of TOM. 

Of the 17k bytes of RAM, Ik bytes were provided on the SBC 80/10A, and 16k 
bytes were provided by an add-on RAM board. Of the 9k bytes of ROM, 8k bytes 
were provided on the SBC 80/lA, and Ik bytes provided on the interface board. 
The memory was allocated as shown in Figure 3.7. 

3.2.2 Interface Hardware — Transferring the analog data into Che 
digital computer and the digital information from the computer to the analog 
outputs requires analog-to-digital converters (ADC) and digital-to-analog 
converters (DAC). An interface board, complete with an ADC and two DACs was 




Baud Rate Factor 
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Parity Disabled 
Parity Gen. = Od 
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Figure 3,5 8251 Mode Configuration 
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Figure 3.4 8251 Conrrand Configuration 
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iustalled in the solar ceil tester. A programmable power supply and its 
digital controller were also installed to allow the microcomputer to step 
through the I-V curve during the measurement period. Each of these devices 
will be discussed in detail. 

The Analog Interface Board The RTI-1200 interface board, raanufac . 'red by 
Analog Devices, Inc., is functionally, electrically, and mechanically 
compatible with the InLei SBC 80/10A microcomputer (3.2). All connections 
between the RTI- ’’00 and SBC 80/10A were made by plugging them into a 
compatible card cage. The RTI-1200 is connected to the microcomputer in a 
configuration commonly called memory mapped interface. The RTI-1200 appears 
to the microcomputer as a block of memory locations. Data and command 
information is transmitted to the interface bcird via software instructions 
that write into memory, and data and status information are received from the 
interface board via software instructions that read from memory. Ths mr^iory 
mapped interface technique simplifies the software instructions needed to 
communicate with the interface board. 

The RTI-1200 performs several functions such as data acquisition, output- 
ting analog signals and memory exp nsion. The RTI-I200’s most basic function 
is data acquisition, as can be seen from Figure 3.8. This is accomplished 
with an analog multiplexer than can accept up to eight differential analog 
input signals, a software programmable gain amplifier, a sample and hol^. 
amplifier, ^nd a 12-bit analog to digital converter. T^e A/D converter can he 
oonfigur-^d for any of three input ranges using wire wrap jumpers. The input 
range of -5V to 4-5V was selected because it pro 'ided the greatest rerolution. 
The programmable gain amplifier was used to increase the input sensit ’ t ivity 
of the vinalog inputs, as shown in Table 3.1. To '•f^ect the gai** netting of 
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Figure 3.8 Analog Interface Board Block Dlagrai 




the programmable gain amplifier, a command word is written to memory location 
address FFF9 hex. The low 2-bits of the command word estabi^ish the gain 
according to Table 3.2. The programmable gain amplifier also enables the 
software program to specify different gains for different channels of input. 

Gain •• 1 Gain • 2 Gain * 4 Gain ■ 8 

-5V to 4-5V -2.5V to '♦•2.5V -1.25V to 1.25V -.615V to .625V 

Table 3.1. Actual Input Range For Various PGA Gains. 

00 - GAIN OF I 

01 - GAIN OF 2 

10 - GAIN OF 4 

11 - GAIN OF 8 

Table 3.2. Command Word to Programmable Gain Amplifier 

Selection of the desired input channel is performed by writing a command word 
to the memory address location FFFA hex. The bottom 3 bits of the command 
word are used to select one of the eight channels, 0 through 7. Once selected 
a channel remains selected until a different channel Is selected. The 
selected channel’s 12-bit input data is read by the microcomputer from memory 
address locations FFFD and FFFE hex. The coding was selected to be two’s 
complement code. 

Outputting analog signals was Implemented using the digital output 
channels provided on the RTI-1200. each analog output channel consist of a 
l2-blt digital to analog converter and an operational amplifier used for 
selecting the output range desired. Digital data Is inputted to the D/A 
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r 

i 

converters by way of memory address locations FFF4 hex through FFF7 hex, as 
I shown in Table 3 >3, using two's complement code. Since the SBC 80/lOA has an 


Memory Byte 

Address Format Name 


FFF7 

( 

0 

0 

0 

0 

SIGN 

BIO 

a 

9 

B 8 

DAC 

1 

LO 

I FFF6 

B 7 

B 6 

B 5 

B 4 

B 3 

B 2 

B 

1 

LSB 

DAC 

1 

HI 

I FFF5 

0 

0 

0 

0 

SIGN 

BIO 

B 

9 

B 8 

DAC 

2 

LO 

1 

FFF4 

B 7 

B 6 

B 5 

B 4 

B 3 

B 2 

B 

1 

LSB 

DAC 

2 

HI 


i 

I 

Table 3.3 D/A Converters Memory Address Locations 

i 

i 

8-bit data bus, the data has to be written to the 0/A converters in two bytes. 
All 12-*bits are simultaneously strobed into the 0/A converter when the four 
I most significant bits are loaded. This is done so chat the 0/A converter's 

analog output will move directly from an old value to the desired new value 

i 

without first going to some intermediate value. 

The RTI-1200 has a PROM socket for memory expansion of 1008 bytes of 
data using the INTEL 2708 Erasable Programmable Read Only Memory (EPROM). The 
bytes are served for the data» control and status functions. 

Programmable Power Supply A programmable power supply was used as an 
! electronic load to obtain the solar cell’s characteristic I-V curve. This was 

j Implemented by Initially setting the current value to zero and Incrementing 

^ the amount of current through the cell while sioiulataneously taking voltage 

1 readings at each of the incremented values, until the entire I-V curve was 

traversed. 

( 

! 

i 
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The ATE 6-5M programmable power supply and the SN 500-121 digital 
progra’^mec, manufactured by Kepco, Inc., were used in the solar cell tester. 
The digital programmer's function was to convert a 12-bit unsigned number from 
the microcomputer into a 0 to 1 volt analog output signal. This 0 to 1 volt 
analog signal was fed into the input of the prograimnable power sup^^y to 
produce a 0 to 5 ampere output current. The computer was then able to control 
the current through the solar cell in increments as small as one milliampere. 

The digital programmer was configured to accept 12-bit parallel data and 
output a single polarity, 0 to 1 volt signal. The parallel input data was 
loaded into a temporary storage register on the negative edge of a strobe 
pulse. The strobe was generated in software, and outputted on port E4H, of 
the microcomputer. A busy signal was then outputted by the digital 
programmer, which remained high until the first clock transition following 
data transfer. Indicating it was ready to accept new data. This signal was 
monitored by input port E5H of the microcomputer. 

The programmable power supply was configured to operate as a current 
source. The output and main feedback capacitors were removed, thereby provid- 
ing the characteristics of a wide band amplifier, ideal for current stabiliza- 
tion and high speed current programming. 

3.2.3 The Solar Simulator -- From a spectral characteristic standpoint 
the most desirable light source with which to measure solar cells is natural 
sunlight, but because of Its variability, solar simulators must be used. 

Three kinds of steady state light sources are commonly used: Xenon arc 

lamps, quartz-halogen tungsten lamps, and ELH lamps (3.3). The Xenon arc 
lamps match the spectrum of natural sunlight closely, but are the most 
expensive of the three and tend to be somewhat unstable. The ELH lamps have 
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an only slightly different spectral response from the Xenon arc lamps, but 
have the advantage of simplicity and lower cost. The quartz-halogen tungsten 
lamps are a very poor choice for they supply too much Infrared light. While 
the ELH lamp Is a quartz-halogen tungsten lamp. It has dlchrolc-coated 
reflectors which transmit a significant portion of the Infrared light while 
reflecting the visible light. The ELH lanq> Is commonly used In 35mm slide 
projectors. There also exists a pulsed type Xenon lamp that has been used 
with automatic test systems to measure array panels. However, pulsed 
simulators are expensive and not particularly stable. A steady state ELH lamp 
simulator, with four ELH lamps, was found quite satisfactory for reliability 
testing. Each of the 300 watt ELH lamps was powered by a 100 vdc power supply 
and cooled with forced air from a fan mounted In the top of the simulator as 
seen In Figure 3.9. A light shutter was used which permitted light to shine 
on the cell only during Che one second measuring Interval, minimizing the 
amount of heat Input to the cell under test so chat It remained essentially at 
ambient temperature during measurement. 

The short circuit current of a solar cell is controlled by the light 
Intensity and Its spectral distribution. While Che ELH solar simulator was 
capable of reproducible results. Its spectral distribution did not perfectly 
match Chat of Che sun. Therefore from Che standpoint of absolute accuracy, 
where ever possible a reference cell of the same type as the cell being meas- 
ured and which had been calibrated In natural sunlight was used to set the 
Intensity. From the standpoint of repeatability It was also Important Chat 
Che simulator's intensity did not drift or fluctuate with time. Individual dc 
regulated power supplies were used for each of the four ELH lamps with a warm 
up period of one hour to minimize drift. A photodiode mounted In the cell 
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holder fixture and maintained at constant temperature was used to nx)nitor the 
intensity between each measurement. The photodiode indicated fluctuation or 
drift in intensity or the presence of a burned out lamp that otherwise would 
have gone unnoticed. 

3.2.4 Cell Holder — In order to position cells under the solar 
simulator in exactly the same place each time and to make stress free voltage 
and current connections, a cell holder was employed (see Figure 3.10). A 
stationary metal ring and an identical shape plastic ring were placed on top 
of each other and clamped together, sandwiching the solar cell leads between 
the rings. This avoided stressing the solder tab when the measuring leads 
were attached. Cells of different size and shape require a different pair of 
rings* A Kelvin four point measuring circuit was used to minimize the effect 
of any external contact resistance* The current contacts were made using the 
vacuum port supports to the back and the metal ring to the top while the 
voltage contacts were made using the thermocouple jacket to the back and a 
clip attached to the solar cell lead. Solar cell characteristics made in this 
fashion will thus show the effect of any contact resistance between the leads 
and the surface metalization and between the surface metalization and the 
silicon, but will be completely insensitive to any contact resistance between 
current carrying elements and the cell leads. Measurements will therefor 
reflect the intrinsic behavior of the cell with its leads, but not variations 
due to connecting the cell to the external circuitry. A foot pedal actuated 
solenoid valve activated the three vacuum ports that held the cell down to the 
thermocouple. 

The cell holder was designed to allow air to be blown over the top and 
bottom surfaces of the cell for temperature control as shown In the photograph 
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Figure 3.9 Diagram of Solar Simulator 
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Figure 3.10 Diagram of Cell Holder 




of Figure 3. LI. An air oilxer box was fabricated which outputs 28° C forced 
air. T .s approach to temperature control has proven to be a big Improvement 
over the previously used heat sink method. Cells with Irregular backs and 
encapsulated cells can be easily maintained at a constant temperature. Figure 
3.12 Is a photograph of the complete hardware system. 

3.3 System Software 

3.3.1 Introduction — Software Instructions for the microcomputer were 
written in a verson of the PL/M programming language, called PL/M-80. This Is 
a high level language written for implementing PL/M on the Intel 8080 micro- 
computer. A high level language was chosen over the lower level assembly 
language because of the ease and speed In programming It offers. PL/M-80 
programming Instructions are oriented more toward the English language than 
toward nemonlcs, resulting In better documentation and a reduction of errors. 
While the PL/M-80 source program can not directly be executed on the 
microcomputer, It can be compiled to a more primitive form called assembly 
language, which can be assembled to create an object program that the 
microcomputer can execute. A program called a compiler exists In the Intel 
Microcomputer Development System at Clemson Unversity and which transforms a 
PL/M-80 source program to an object program. The development system also 
provides a program culled the Intel Systems Implementation Supervisor, ISIS-II 
which gives the user the capabilities to do text editing, linking, and program 
locating, thus reducing the software development and debugging time. 

With PL/M-80, programs should be writ tin in modules. Each module can 
then be checked Independently for errors before being linked together. This 
localizes errors for easier debugging. The software for the tester has one 
Main Program Module and a number of other modules containing only procedures. 
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Figure 3.11 Photograph of Ceil Holder 
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Each ot these procedures and the Main Program Module will be described In 
detail. Detailed flow charts and the complete PL/M-SO source program are 
contained In Appendix B. 

3.3.2 Main Program Module The Main Program performs a sequence of 
call Instructions. Initially the cell Identification procedure, named the ID 
procedure, Is called to Instruct the operator to enter the cell type, lot 
numer, aud the stress level. The Main Program remains In a loop with each 
pass through the loop being a measuriment, display, and output of the para** 
meter data for one solar cell. 

3.3.3 ID Procedure — This procedure Is called from the Main Program 
whenever a change In cell lot Is needed. The operator Is Instructed to type 
the cell type, lot number, and stress level. This Information Is stored to be 
later punched on paper tape along with the parametric data of the cell. 

Before the procedure returns to the Main Program It types out the table 
headings for the parameters. 

3.3.4 Get Cell Number Procedure — This procedure Is called by the Main 
Program. Initially the procedure prompts the operator to type in the cell 
number, by printing 'CELL NUMBER ■ '. The procedure will accept up to three 
characters and store them in consecutive otemory locations. They will be 
printed out when the cells' parameters are printed. In order to correct any 
character that was pressed In error, only the lest three numbers are saved. 
After entering the cell number a carriage return will return the procedure 
back to the main program. 

3.3.5 Step Through Procedure — This procedure is called from the Main 
Program approximately one half second after the light shutter is activated to 
open. This delay, which can be varied, allows the light shutter to be fully 
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open before any data points are measured. It then stepw through the 
characteristic I-V curve by repetitively callng the Ccnipute Current Procedure 
which outputs the next value for current, and upon returning measures the 
corresponding voltage across the cell. Thi. process Is repeated until the 
cell has passed completely through the I-V curve. Since the programmable 
power supply has a finite settling time a delay of about ten milliseconds 
occurs between a change In output current and the time Che voltage Is read. 

The voltage Is then read repetitively until two consecutive voltage readings 
are the same. Since Che same voltage potential could be measured on the rise 
and fall of an oscillation the initial 10 ms time delay was included as an 
extra precaution. Both the delay and voltage comparison readings are used to 
ensure the current has settled down. If the voltage readings have not settled 
down after sixteen times the program stores the last value and continues to 
the next data point. This situation may occur In the flat region of the I-V 
curve, close to the short circuit current, where a small amount of ripple on 
Che programmable power supply's output current can cause a larger ripple on 
Che voltage. The inability to have the voltage settle down in this region of 
the I-V curve does not significantly affect measurement of the short circuit 
current parameter. The procedure returns to Che main program upon a negative 
voltage reading, which Indicates the solar cell has entered the second 
quadrant of the I-V curve. 

3.3.6 Compute Current Procedure — This procedure Is called from the 
Step Through Procedure. It determines the value for current of t*-e next data 
point along the I-V curve, ouputs It uO the programmable power supply and 
returns to the Step Through Procedure where the corresponding voltage value is 
read. 
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There are three regions of the I-V curve which require special concern 
to ensure that a sufficient number of data points are obtained in order to 
accurately determine cell parameters* Data points measured in the vicinity of 
the open circuit voltage are used to calculte the slope dV/dl at I » 0 ^hich 
is proportional to the series resistance^ and to determine the open 
circuit voltage, Data points measured in the vicinity of the maximum 

power point need to be closely spaced in order to calculate the maximum power 
point, P^, accurately. Data points in the vicinity of the short circuit 
current, Ig^, are used to determine both Ig^ and the slope dV/dl at 
V * 0 which is proportional to the shunt resistance, Therefore the 

Compute Current Procedure has separate modes of operation for each of the 
three regions in order to ensure that sufficient data points are measured to 
decermine the cell parameters. Figure 3.13 illustrates the regions in which 
each mode is applicable 

The Compute Current Procedure is initially in mode 0, with an output 
current of zero asperes. In mode 0 the current is incremented by 32 milll- 
amperes each time the procedure is called, resulting in another data point. 
Each time a data point is measured, the power output is calculated by 
computing the product of current and voltage. This value for power is 
compared with the value from the previous data point to determine if the 
maximum power point of the cell has been reached. If the power from the 
previous data point is greater, then the maximum power point has been passed. 
The procedure then goes into mode 1 where the current is decremented by 1 
milllarap for each data point. The power of each data point is again 
calculated and compared to the previous data point. Mode 1 continues until 
the calculated power decreases indicating that the maxiniura power point has 
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again been passed* The previous data point is then identified and saved as 
the maximum power point and the compute current procedure changes to mode 2* 
Mode 2 initially outputs the current value that mode 0 ended with followed by 
a value incremented by 4 ma. The change in voltage between these two data 
points is used to determine the amount of increment to be added to the current 
for the next point. The change in voltage is again calculated and the next 
increment determined and so on for the remainder of the curve. This procedure 
is necessary rather than using a fixed increment because of the large range of 
possible slopes in the region around the short circuit current. Not only does 
the slope vary from cell type to cell type, but it also varies within cell 
types as a result of degradation due to stress. When measuring a high 
efficiency unstressed cell with a flat slope, dl/dV near zero, an increment as 
small as one milliampere is needed to obtain a sufficient number of data 
points in the vicinity of Ig^, where as a badly degraded cell requires 
larger increments for the same number of data points. Use of the 1 ma 
increment on the degraded cell would result in an excessively large number of 
data points. Mode 2 uses the voltage change between the two data points to 
determine the size of the increment. If the change in voltage for the 4 ma 
increment is between 1.2 mv and 19.5 mv then the increment remains at 4 ma. 

If the change in voltage is larger than 19.5 mv then the increment is reduced 
by 1 ma while if the change in voltage is below 1.2 mv the increment is 
increased by 1 ma. (Ima is the minimum increment possible). The Compute 
Current Procedure remains in mode 2 until the I-V curve has been completely 
measured (voltage becomes negative) and is reset to mode 0 for the next cell 
to be measured. 
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Figure 3.13 Measured Cell Characteristic Illustrating 
Different Compute Current Mode Regions 
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3 . 3-7 Scope Procedure — The Scope Procedure is called from the Main 


Program after all I*V data points have been measured and stored in a table in 
memory* The scope procedure executes in a loop that outputs the voltage and 
current values of each data point to D/A converter one and D/A converter two 
respectively* The outputs of the D/A converters are connected to the X and Y 
inputs of an oscilloscope to display the I-V curve. Return to the Main 
Program occurs after any key is pressed on the TTY keyboard. 

3.3.8 Print Procedure — This procedure outputs a string of characters 
to the TTY. The characters need to be in consecutive memory address 
locations, and in their ASCII representation. The memory location of the 
string of characters and the number of characters to be printed are the 
parameters that need to be passed to the print procedure when it is called. 

3.3.9 Fetch Procedure — The Fetch Procedure inputs a string of 
characters from the TTY and stores them in consecutive memory address 
locations. The starting memory address location and the number of characters 
to be input need to be passed to the procedure when it is called* 

3.3.10 Numout Procedure — The numout procedure converts a number from 
binary form to an ASCII string suitable for printing. This procedure requires 
the following parameters when called. 

VALUE The number whose printed representation is desired. 

BASE A n integer from 2 to 16 specifing the base of ’value*. 

LC Leading character, either a * ’O’, or * 0 *. 

BU^AD R Buffer address. 

WIDTH T he number of characters desired in the printed 

representation. 
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3.3.11 Tbout Procedure ~ This procedure outputs 12-bits of data to the 


programmable power supply's digital controller. The lower 8-bits are 
outputted on port 2A, and the upper 4-bits are outputted on port 2B. A strobe 
pulse outputted on port lA latches the 12-bits of data into the controller. A 
busy line received from the controller is monitored on input port 2C until it 
goes low, indicating the controller is ready to accept more data. The 
procedure returns to the program from which it was called when the busy line 
goes low. 

3.3.12 X1232 Procdure — This procedure takes a 12-bt signed Integer 
and converts it to a 32-bit signed integer in order to convert the voltage and 
current values read from the A/D converters to a format acceptable for the 
floating point arithmetic library routines. It can be called from any proce- 
dure that declares it external, provided the memory location address of the 
12-bit number to be converted and the low memory location address where the 
32-bit number is to be stored are given. 

3.3.13 Pmax Procedure — This procedure is called from the Step Through 
Procedure. It calculates the product of the voltage and current values of the 
last data point measured. This value for power is compared to the previous 
data point’s power output. A flag, LT, is set to one if the ^ast data point’s 
power becomes less than the previous point’s, otherwise the flag remains equal 
to zero. Before ^ Ae procedure returns to the Step Through Procedure, the 
largest power point value is stored in memory. 

3.3.14 Global Variables Procedure — This procedure was used to declare 
public variables that were used in more than one procedure. The variables 
were declared external in each of the other procedures. This was used to 
simplify the bookkeeping so as to reduce the possibility of declaring a varia- 
ble public in more than one procedure. 
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3.4 Operating Procedure 


3.4.1 Power Up — In order to prepare the system for mo.asurement the 
following power up procedure Is used: 

1. Turn the main power switch to the ON position. 

2* One at a time, set each of the four £LH lamp switches to the ON 
position. There is a built in delay of approximately 3 seconds and the opera- 
tor should note that each lamp is operating before turning the next on. T’le 
lamps should be allowed to operate for 1 hour prior to making measurements. 

3. Turn the teletype power switch co the LINE position. 

3.4.2 Start Up — In order to begin measurement taking the following 
start up procedure is used: 

1. Press the RESET button located on the front panel. This resets the 

microcomputer and initiates execution of the system monitor, a program 

used for debugging and program start up. A message, **BLC SOP MONITOR REV.C** 
is printed out Indicating the system is ready for a monitor command. 

2. To start execution of the main tester program the user types: 

”G580 )*’ ( ) Indicates a carriage return). A is printed out indicating 
that the main tester program is waiting for a command. 

3. One of two valid commands may now be entered: a **0 y* to start the 

routine IV characteristic measurement program or a "1 y to start the time 
varying effect measurement program. If an invalid command is given a ”?** will 
be typed out and the correct command should be entered. 

3.4. Routine IV Measurement Program — Having entered ’*0 the user is 
prompted to type in the cell type, cell lot number, and stress level as 
follows: 
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CELL TYPE - B ^ 

LOT # « 19 > 

STRESS LEVEL - 3 > 

where underline indicates information entered by user. Table headings are 
then typed out: 

TYPE*LOT*S * VOC * Ig^ * * '^m * ^m * TEMP* 

The user is then prompted to type in the cell number and temperature as 
follows: 

CELL NUM - 35 > 

CELL TEMP - 28,0 ^ 

Ultimately it is planned to receive the temperature Information directly from 
the thermocouple readout, but at the moment it must be manually entered* 

After pressing the carriage return upon entering the temperature data, the 
pneumatically controlled light shutter opens, the IV characteristic is 
measured, and the shutter closes. The sequence takes less than 1 second. The 
digitized IV curve is automatically displayed on the monitor oscilloscope. If 
the curve appears satisfactory, pressing any TTY key will cause the measured 
parametric data to be typed out in the previously prepared table and 
simultaneously punch the informaton on paper tape. If the IV curve is not 
satisfactory, holding down the OFF button on the paper tape while pressing any 
TTY key will erase the data without recording. 

After the parametric data has been printed out and the paper tape punched 
the prog’^am will have executed one pass through a loop and returned to prompt 
the user to type in the next cell number. The program will remain in this 
loop until all the cells in a lot have been measured. The operator can then 
press the RESET button and begin execution of the main tester program again. 
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3.4.4 Time Varying Effect Measurement Program — Having entered "1 >" 


the user Is prompted to type In the number of measurement!^ cell type, cell 
lot number, and stress level as follows: 

# OF MEAS - 20 ^ 

CELL TYPE -li 
LOT # - 12 > 

STRESS LEVEL - 3 > 

After the last carriage, table headings are typed out as before 
TYPE * LOT * S * VOC * Isc * Pm * Vm * Im * TEMP* 

The user Is then prompted to type In the cell number and cell temperature as 
follows: 

CELL NUM - 25 > 

CELL TEMP - 28.0 > 

After pressing the final carriage return the pneumatic light shutter opens, 
the IV characteristic Is measured and the digitized data displayed on the 
monitor. At this point the shutter does not close, but remains open. Another 
measurement will then be made, stored In memory, and displayed on the monitor 
every time a TTY key Is pressed. Usually It Is desired to make measurements 
at regular Intervals. After the final measurement the parametric data Is 
printed out In the previously prepared table and curves plotted. To ic asure 
another cell the RESET button Is pushed and the process started again. 

The time varying effect measurement program makes It possible to measure 
changes in parameters that occur as a function of illumination time. Specific 
examples will be discussed in the next sectlor . It should be noted, however, 
that with the shutter open the temperature will rise and asymptotically 
approach a higher value. In this mode of operation it is necessary to record 
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Che temperature of each run and to use this value to correct Che data. 
Presently this must be done manually, but ultimately when tempe;:ature data is 
directly accessed this will be automatic. Thus at present the temperature 
Chat is entered in Che prompting mode and Chat is displayed in Che data table 
is the starting (dark) temperature. 

3.5 Time Varying Effects 

Use of the short interval tester uncovered a parametric variation with 

time of illumination, which had gone unnoticed using conventinal test methods. 

Figures 3.14 and 3.15 show typical variations of V , I , and P with lllu- 

oc sc m 

minaclon time. This type of change has been reported previously as photon 
degredation (2.4). This time varying effect was observed on every cell of one 
n/p type cell and on a few cells of another n/p type. All cells recovered to 
their original performance after being in the dark for a few minutes. 

Figure 3.16 shows plots from the short interval tester of one susceptible 
cell's output Initially, after an illumination of 20 seconds, and after 10 
minutes. The tester can be a useful tool for further study of this 
phenomenon. Digitized I-V curves can be measured as often as every second and 
scored in memory to be plotted and/or printed out later. Reduced parametric 
data can also be printed out. 

While the time varying phenomenon is interesting and deserves further 
study, Che variability of Che cells pose some problems when crying to measure 
these cells accurately. Because of the non-repeatability when making measure- 
ments' in close succession, the standard cell used for adjusting the light 
to one sun illumination had to be chosen from a non-time varying cell type. 
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F\gure 3.14 V and I as a Function of Illuminations 
Time for Representative Cells 


111 


Figure 3.15 as a Function of I 
for Representative Ce 
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Figure 3.16 Voltage-Current Plots for Different 
Illumination Times 
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rhe question also arose of whether to p. : Iluniinate the cell and wait until 
thr ’-^gradation occurred before measurin;*;. to measure immediately after 
i i iuminatioHM It was found that differc^wr.s as much as 5% in 4% in 

and 9 % in Pjjj r3sulted between th . L*'ethods of measurement. 

Measu/*nf^ immediately upon lllumi ..as chosen because repeatability was 
of greater concern than abs'^lur.. .. 4, acy although the actual power tended to 

be overestimated in this appro.u Also, a geat deal of time and effort would 

be required to pre-illumlnate rhe cells before measurement. 

3.6 Data Management 

In order to statistically analyze the electrical data obtained on cells 
subjected to accelerated stressing it is necessary to enter the parameter 
values, cell number and type, and the stress conditions into an IBM 370 
computer. The data management routine which accomplished this prior to the 
short interval tester, is shown in Figure 3.17. The output of the original 
test procedure was hard copy of the IV characteristic curve and values copied 
from meter readings of 1^^, and temperature. It was necessary to first 

obtain the maximum power point from the characteristic using an overlay of 
constant power hyperbolas. The coordinates of the maximum power point 
together with all the other measurement parameters were then transcribed onto 
keypunch forms. Cards were then prepared from the information on the form and 
the dat;' placed on disk for permanent storage and later manipulation. Errors 
were introduced in graphically obtaining the maximum power point, in 
transcribing the data, and in keypunching the data. Construction of the short 
interval tester provided an opportunity to eliminate these errors. 




DATA ANALYSIS 


Figure J.17. Initial Method of Transferring Measurement 
Data to the IBM- 370 Computer 
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The present data management procedure used with the short interval tester 

is shown in Figure 3.18. The tester stores the voltage-current data pairs 

taken on each cell in random access memory. It also stores specific values of 

V , I , I , V , and P as well as the manually entered test conditions. The 
oc sc m m m 

specific parametric values are then recorded on paper tape, but because of the 
bulk of this storage medium it is not practical to save actual VI data point 
coordinates and these are discarded. The paper tape is then processed on a 
PDP-15 where the parametric values are transcribed onto magnetic tape in an 
IBM compatible format- Use of the PDP-15 is necessary since the IBM equipment 
does not have paper tape reading capability. The PDP-15 generated magnetic 
tape can then be manually transported to the central computer facility where 
it is read onto an IBM-370 magnetic tape for permanent storage. Whenever data 
manipulation is required the information is placed on disk. Permanent storage 
of data on tape instead of disk results in appreciable cost savings. This 
particular data management system eliminates any direct human intervention, 
but is somewhat awkward in that several physical transport steps are involved 
between different computers. 

Introduction of an IBM compatible floppy disk into the short interval 
tester system will eliminate the need for paper tape as a storage medium and 
consequently the need for the PDP-15 computer conversion to magnetic tape as 
shown by the dashed line path of Figure 3.18. Up to 300K bytes of data can be 
stored on a single sided diskette. Not only could they be used to store cell 
parametric values, but they could also store all the actual current-voltage 
data points so that the actual V-I characteristic curve could be reproduced 
later. In addition, software could also be stored on diskettes making it 
easier to update the system and add programs. Since almost all the PROM 
memory area is presently used program additions will require either an 
additional PROM board or a floppy disk system. 
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SHORT INTERVAL TESTER MEMORY (RAM) 



Figure 3.18. Present Method of Transferring Short Interval 
Tester Data to IBM-370 Computer 
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If the need arises in the future for a substantial number of hard copies 
of IV characteristics, a graphics system can be incorporated into the system 
replacing both the oscilloscope and the x-y plotter* 

3.7 Cost 

As microcomputers drop in price and improve in performance they are being 
used in applications where the larger minicomputers were used previously. In 
the case of the short interval tester a system was needed to perform dedicated 
functions and the lower cost microcomputer was ideal* Since the system was 
built in house, the total cost was reasonably low. Table 3.4 shows the cost 
of the major capital equipment items required to construct the system. 


Microcomputer 

$ 

1600 

Analog Interface Board 

$ 

1000 

RAM Memory Board (16k x 8) 

$ 

300 

Programmable Power Supply 

$ 

600 

Digital Controller 

$ 

400 

Teletype 

$ 

1400 

Plotter 

$ 

1700 

Oscilloscope 

$ 

800 


$ 8300 

Table 3.4. Cost of Major Capital Equipment 

A major consideration in the design of the system was minimizing capital 
equipment cost. For this reason several peripheral such as the Teletype and 
X-Y plotter were used instead of a more satisfactory CRT, line printer and 
digltvil plotter combination. Also a floppy disk system was omitted in the 
original design. These Items could be added to the system later as money 
becomes available to Improve the system*s capabilities. 
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Another major cost to a computer system is the cost of software. While 
the cost of the hardware has decreased due to technolgy the c'^st of the 
software has risen as labor rates increase. The cost of software can be a 
significant portion of the total first cost. In part the Intel 809/10 
microcomputer was chosen because of the facilities at Clerason for writing 
software in a high level language, PLM/80, which saved both writing and 
debugging time. Even so, it is estimated that to reproduce the software that 
has been developed thus far for the tester would require in excess of $100,000 
in an industrial environment. 
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4.0 ANTIREFLECTIVE COATING EXPERIMENTATION — H.A. Walker 


The possibility that the electrical characteristics of photovoltaic 

cells degrade as the intiref lective (AR) coating is removed was theorized 

during earlier testing cf F, G, and H cells. (See Second Annual Report, April 

1980.) An an attempt to quantify this, new visual inspections were made of 

each cell which had undergone 500 hours of pressure cooker stress and an 

estimation made of the percentage of AR coating that had been leached away. 

These estimates were compared with the percent decrease of both I and P . 

sc m 

Figures 4. 1-4. 3 show the estimated AR coating missing vs. thepercent decrease 

of I and P for each single cell. The plots do not all show a definite 
sc m 

trend, but the G cells, which contained the least AR coating, seemed to 
exhibit a greater degree of correlation than either the F or H cells. Also, G 
cells are identical to I cells except for the absence of external leads. 
Therefore, these findings, though baseu on rough estimates, were significant 
enough to justify further research. 

Mr. Edward Murphy of MIT Lincoln Labs first suggested the possibility of 
quantifying AR coating loss through colorspectrum analysis. We were fortunate 
in having available an IBM 7400 spectrophotometer which could be used for 
this purpose. The 7400 is a fast-scanning, high-precision, stable 
spectophotoraeter that operates in the visible wavelength range of 400 to 70 
nonometers, and measures the light reflected from a sample at designated 
wavelength intervals (4.1). The light is filtered through a rotating circular 
interference filter assembly to provide continuous wavelength selection. 
Reflected light is fed to the color sensor by fiber optic cables from the 
diffuse reflectance sample holder. The sample holder contains a fixed fiber 
optic sensing head and a sample pedestal which can be adjusted vertically. 
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Visuslly EscitnaCfid • of AR Coacintt Rctnovad vs. Measured 
% Decrease in Isc •’m Lots F-14 and F-19 after 
500 Hours Pressure Cooker 
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PERCENT DECREASE IN ELECTRICAL OUTPUT 


Figure a . 3 Visually E.^riTnaCed % of AR Coating Removed vs. Measured 

Decrease in Tsc and Pm for Lots H-14 and H-19 after 
500 Hours Pressure Cooker 
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The area measured is a circle with a diameter of 32mm. A solar cell could be 


placed in the pedestal and locked into position directly below the sensing 
head. 

Control of the instrument is accomplished through a Color Sensor 
Analyzer which also computes reflectance data and provides hard copy color 
information. The spectral relfectance data obtained from the 7400 system is 
used to calculate color space coordinates, color difference, chroraasticlty 
coordinate values, and cristimulus values. (See Appendix C for a discussion 
of these terms.) 

To test the feasibility of using this instrumentation another trial was 
held using a small group of C cells which had undergone various color changes 
as a result of pressure cooker stress. On some cells onlu a change in hue was 
observed, while others exhibited a substantial loss of AR coating. The 
resulting wide variety of colors suggested a good test of the IBM system’s 
apparent visual accuracy. Various locations on each cell were analyzed and 
the tristimulus values and color space coordinates measured. Favorable 
comparisons were observed between these readings and the general visual 
appearance of each cell to the eye. Some concern had been expressed that the 
surface metallization would have a considerable effect on the measurements- 
Tc test this possibility, readings were taken on two cells before and after 
painting their metallization black. The only significant effect noted was a 
decrease In the overall lightness (L color space coordinate) of the cells, 
along with equal percentage drops in each tristimulus parameter. 

The next step was to measure untested cells and to follow them through 
the stress testing procedure, measuring them at each down time. I-cells were 
used and three readings were taken for each cell. Since the I cells had not 


been stressed, their readings were more uniform over a given cell. 

Measurement repeatability was the main concern, so several sets of 
measurements were made over a period of several days to determine the 
Importance of exact cell positioning. These tests determined that it was 
necessary not only to measure the Identical area each time, but also to place 
the cell on the pedestal in exactly the same orientation each time. This was 
accomplished by fitting an allignment mask over the cell as shown in Figure 
4.4 so that only the desired area was seen by the measuring sensor. 

Orientation of the cell to the pedestal was maintained constant by lining up 
the center solder line on the cell with the edge of the sample holder. A 
similar mask was placed over the cell during electrical measurements so that 
only the areas for which the color readings were taken, would affect the 
output. These methods resulted in consecutive color measurements (of the same 
cell area on dilfferent days) that agreed within two percent. 

Another factor to be considered was observer angle. The 1931 CIE stands 
ard angle referred to a 2 ° visual field and dealt mainly with distant viewing, 
whereas the 1964, 10® standard was set for relatively close observations. To 
investigate possible differences between the standards, two sets of readings 
were made for each cell. Only slight differences were noted between measure- 
ments of the same area. Therefore, the 10® angle was used for all susequent 
spectral analysis. 

The measured I cells are currently undergoing pressure cooker tests and 
results will be reported in the next annual report. 
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Figure 4.4 Alignment Mask in Position on an I-Cell. 
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CONCLUSIONS 


5 • 1 D evelopment of a Mathematical Model to Assess the Reliability 
Consequences of Operation in the Second Quadrant 

A computer model based on the heat flow properties of silicon and the 
encapsulation materials, and on the variation of the electrical properties of 
silicon as a function of temperature, has been developed which relates ceil 
tenperature to cell VI characteristics in the second quadrant- This model 
predicts that as a cell’s reverse bias voltag- increased the cell 
temperature will increase, uniformly at first, and then non-unif ormly with 
suosequent hot spot formation as the voltage approaches breakdown* Agreement 
with experiment is excellent for unencapsulated cells and reasonable for 
encapsulated cells* Model parameters are less well known for encapsulated 
cells and the measurement problems greater, which accounts for the differences 
observed with encapsulated cells* Nevertheless, the model provides the only 
reasonable method of estimating cell temperatures in modules where cells 
become reverse biased* There seemed little difference between encapsulatea 
and unencapsulated cells regarding the values of S and f , the fitting 
parameters, althrough the sample size consisted of only 8 cells, 4 
encapsulated and 4 unencapsulated. For most cells S'* 0*9 + 0.2 and y ■ 0*3 + 
0*1* Additional cell types should be investigated to see if Sandy vary from 
structure to structure. 

The ability to relate cell temperature to electrical characteristics can 
provide the module designer with valuable guidelines. For example, experi- 
ments with EVA encapsulation have indicated it would not withstand tempera 
tures in excess of 140^C* The model predicts this temperature will be reached 
in glass superstrata modules at 15 volts and in masonite modules at 10 volts 
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and therefore under no circumstances should the voltage be allowed to exceed 
these values. The other components of the encapsulant and interconnection 
systems such as the glass, masor*ite and solder also have their own critical 
conditions. For glass the time rate of change of temperature may be the 
critical parameter, for masonite, the evolution of moisture, and for solder, 
its melting temperature. Rather than examine each component at this point a 
general statement may be made that cell temperatures should be kept as low as 
possible by restricting voltage. Based on the results of both our limited 
miniraodule work and our extensive unencapsulated cell accelerated stress 
program this means cell temperatures should be kept below 100®C, which for 
most cell types translates to voltages less than 7 volts. This would mean 
bypass diodes for at least ev^ery 14 cells. Jlore frequent bypass diodes are 
preferrable aid the new technology section describes a design which was formu- 
lated to fabricate Integral bypass diodes on every cell. 

5.2 Development of Improved Solar Cell Test Methods 

The microcomputer controlled, constant tempercture, short interval 
tester described in this report is a great improvement over the previously 
used manual method of plotting the I-V curve, reading and recording data 
points- It has decreased both the measurement and data reduction time by a 
factor of 5 or more, along with essentially eliminating human error. At the 
same time the user retains a feel for the measurements taking place because of 
the monitor oscilloscope that displays the measured data points. 

Accuracy and repeatibillty have been improved by measuring cells at a 
more nearly constant temperature. The fraction of a second required to obtain 
the data enables the light source to shine on the cell for such a short time 



that the heat generated by the light source has little affect on the cell's 
temperature. 

Time varying effects may now be observed and quantified using the new 
ueasureraent system. It is not clear how stressing affects a cell's 
sensitivity to the phenomenon, but this can now be investigated. 

5.3 Color Measurement of AR Coating 

Commercial color measurement instrumentation can be used to characterize 
AR coatings reproduc Ibly . While this characteriz tion is appropriate only to 
the visible portion of the spectrum evidence suggests that color may be rela- 
ted to the cell's PV output. A series of quantitative tests is underway to 
determine if this is the case. 
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6.0 RECOMMENDATIONS 


6.1 Improved Test Methods 

The capability to directly access temperature data should be incorpor* 
ated in the tester. A floppy disk should be added to simplify the transfer of 
data to the IBM 370 system. A floppy disk would greatly increase the general 
versatility of the system and permit additional tasks to be performed. 

6.2 Second Quandrant Model 

Improved temperature measurements as a function of reverse bias need to 
be made on a variety of different cell types under different illumination con- 
ditions , both encapsulated and unencapsulated. This will enable the fitting 
parameter used in the model to be better defined which in turn will allow an 
assessment of the temperature rise anticipated in field installations under 
reverse bias conditions. Together with accelerated test data this will permit 
an estimate of array power degradation over life to be made. 
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7,0 NEW TECHNOLOGY 


Figure 7.1 illustrates a solar cell with an integrated design of bypass 
diode. This design features a wrap-around top contact, a back contact in the 
center of the cell and an integrated bypass diode. The heavy line in the 
schematic in Figure 7.1 represents the preformed V-shaped interconnect in the 
isometric view. Since it only connects to ore side of the cell this metal 
lead could easily be attached by automated equipment. 

It is anticipated that the diode and the cell junction would be manufac- 
tured in the same diffusion step. The diffused layer would then be removed 
from around the diode and from the back contact area. After this metalization 
would be applied and the interconnect preform soldered using essentially the 
same technology as for state of the art cells. 

In order to illustrate the advantage of a one diode per cell 
configuration a module of 36 cells connected in series was fabricated. This 
module was constructed so that it could be operated either with or without 
bypass diodes which were fabricated using solar cell material and processing 
to simulate the characterisitic expected from integrated diodes. Figure 7.2 
shows the VI characteristics of the module under different shading conditions. 
While the exact amount of power loss experienced by the module will depend on 
the module load it is clear that the loss will be much greater without diodes 
than with. 
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Figure 7.1 Proposal Solar Cell with Integrated Bypass Diode and Wrap Around Top Contact 
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Figure 7.2 Output Characterise ic of 3b Cell Module With 
and Wltout Bypass Diodes. 
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APPENDIX A. MATHEMATICAL TOOLS 


The thermal balance equations derived in Section 3.3 are of the form 

2 


_1 3_u 

3x2 X 


k(u) du 9x ^ ^ k(u) 9c* ^ 


and are non-linear, parabolic partial differential equations. An excellent 
treatise on these and other similar types of differential equations is written 
by D.U. von Rcsenburg (3.16), **Methods for the Numerical Solution of Partial 
Differential Equations". 

The Crank-Nlcholson technique was the method used for evaluation. In 
this technique all the finite differences are written about the point, x^, 

half way between the known and the unknown time loads. In Figure A-1 
this point is shown as a cross while the black dots denote known values of 
the dependent variable. Values of the dependent variable, u, are computed 
only at the points designated by open circles. The second-order-correct 
analog of the time deri^^ative at the point x^, is; 


jji 

V t i , n+4 
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i,n 
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A geometric interpretation of the analog is shown in Figure A-2. Thus the 
time derivative is approximately given by, 


u. u. 

3 ji 1 , n-H - I , n 

3 c i , n+^2 t 


(A-3) 


2 

The approximation of the second derivative, , witnout requiring the 
evaluation of the dependent variable at the time level for which the time 
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DEPENDENT VARIi 


Figure \< 



TIME 


!. Geometric Int erpre tat ion of a Socond-Order-Cor rec t 
First Derivative 


index has the value of n+^ Is the key to the Crank-Nicholson technique. This 
derivative is approximated by the artihmetic average of the finite difference 
analogs at the points x^ , t^ and , ^n+ 1 ' 


,n "^i ^n 

. . I 2 , ,2 

i,n+S (^x) 


The first derivative is approximated 
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3x 


i , n+S 


^ 1 ^i-l,n 

2 2Ax 


i ^n ^ ,n^I 

in the same manner: 
2 Ax 


^ ^1- 1 ,n-H . ( 

(Ax)* 


(A-5) 


This leaves the evaluation of the coefficients of the equation at the point 

X. > t and the nroblems with the term ~ . The coefficients can be 

1 n-fU ‘ 3x 

determined by an iteration procedure in which the time dependent value is first 

calculated using the value of the coefficient at the point x., t . This gives 

i n 

a value for u at x^, ^'n+1* which the equation is solved again by using 

values for the coefficients that are calculated by the use of the arithmetic 
average for u at x., t and x., t ... Next the time dependent value is 
calculated again, and so on, until a predetermined accuracy is reached. This 
method was found to be unnecessary, however. Instead, the value at x^, was 
used and the time step was taken small to make up for the inherent danger of 
oscillation. In the same manner the square of the first space derivative was 
handled. If equation (A-5) were used it would have been necessary to solve a 
matrix of se ond order equations, which would have been very time consuming, 
if not impossible. 

Using the described method the following set of o<piations was obtained; 
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A set of equations of this sort, a tridiagonal matrix, is readily solved by a 
computer with the algorithm of Thomas. The algorithm Is as follows: 

first compute, S. * b. r with & = b , 

^ I o o 
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to obtain the dependent variable with 
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n n+l 
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and 




(* from previous period.) 

The solution of equation (A-1) can be approximated using the methods 
described, but care must be taken to limit the time step in order to suppress 
instabilities caused by the coefficients and by the square of the first 
spatial derivative. It is not necessary to take an excessive number of 
intervals along the spatial axis however. 
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The core of the program is the subroutine that calculates the 
temperatures at the point t. ^ . Figure A-3 lists the flow chart of this 
subroutine for the unencapsulated cell. It is largely self explanatory except 
for the boundary conditions which are discussed In Section 2.5. The subroutine 
contains two limiters: a temperature limiter and a current limiter. The 

temperature limiter is necessary to reduce the risk of overflow, as for 
example when the exponential functions in the current density equation are 
calculated. The current limiter is essential for physical reasons. As 
mentioned earlier, a solar cell is essentially a current source, and its 
operation in the second quadrant will typically be current limited rather than 
voltage limited. A secondary reason for the use of a current limiter is 
computational; the I^V curve of a cell in the second quadrant is single valued 
in voltage , but often double valued in current. 

The computational current limiter works in a similar fashion to a 
physical current limiter. If current exceeds a preset value, the voltage is 
reduced until the current falls to the present value. The current limiter has 
a typical overshoot of about 5 computational periods. Both the temperature 
and current limiters set a flag when they are used. 

Figure A-4 shows the structure of the complete program. As mentioned, 
the program solves the I-V curve by stepping along the current axis. For each 
current level the program solves the steady state solution with decreasl ig 
time interval for fast converquence . The program is capable in solving 
transient solutions, but should use a fixed, small time interval to be 
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NEW TEMPERATURE CALCULATION 


I A= 1 
1=0 

TJ=TB( lA) 

CALCULATION OF CURREN TDCNS 1 I Y 
CALL CUROEN 

CALCULATION OF THERMAL CONDUCTION 
WK=WSI/( .58+TJ/350) 

CALCULATION OF THERMAL DIFFUSION 
A=Wk/( 1.6H7*WSI ) 

CALCULATION OF FIRST SET OF MATRIX ELEMENTS 
Al ( IA)=0 

Bl( IA) = -2/0R**2-1/( A*DT) 

Cl ( IA)=2/DR«*2 

OK IA) = (2/DR**2-T/(A*0T ) )*TJ-2*TB( I AM )/DR**2 
C- ( C0*V-HL«TJ+HL**T I NF+ I RR )/WK 
WRITE(6.10)AI(IA),BI(IA),CI(IA).DI(IA) 

00 300 1=1,74,1 

IA= 1+1 
1J=TB( lA) 

CALCULATION OF CURRENTDLNS I lY 
CALL CURUEN 

CALCULATION OF THERMAL CONDUCTION 
Wk=WSl/( .58+TJ/350) 

CALCULATION OF THERMAL DIFFUSION 
A=WK/( 1 ,647«WSI ) 

CALCULATION OF THE REMAINDER OF THE MATRIX ELEMENTS 

Al ( IA) = ( 1-1 )/(( 2*1-1 )*DR**2 ) 

HI ( IA)=-l/0R**2-l/( A*Or ) 

Cl ( I A| = I / ((,"'* I -1 !*0R**2 1 

01 ( IA) = -AI( IA)*TB( IA-l) + ( l/DR**2-l/( A*DT ) )*IJ 
C-C I ( I A 1*1H( I A+n-( CO*V-Hl •TJ+Hl *T I Nl + I RR ) /Wk 
C*( UT( IA+1 )-1B( IA-1 ) )**2/( (205+TJ )*4*DK**cM 
C WR Ml ( 6 . 1 0 ) A K I A ) . B I ( I A ) , C I ( I A ) , t) I ( I A ) 

300 CONMNUE 


Figure A-3. Temperature Cdlculations Subroutine 
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I - lb 
I A= I '►I 
TJ-- U3( I A) 

CALCULATION OF CURRENTDLNS I TY 
CALL CURDEN 

CALCULATION OF THERMAL CONDUCTION 
Wk=WSI/( .58-^^TJ/350) 

CALCULATION OF THERMAL DIFFUSION 
A=Wk/( 1 .6U7^WSI ) 

CALCULATION OF THE 76TH MATRIX ELEMENTS 

Al ( IA) = 1/(2^0R^^?) 

8l( IA)=-1/(2^0R»*2)-l/( A^DT ) 

Cl ( I A)=0.0 

Dl ( IA)=-AI ( IA)*TB( lA-1 )-(BI{ IA)+2/(A*0T) )^TJ 
C-( CD^V-HL^^TJ-mL'^T INF+ IRR)/WK 

C+7S. S^^HL^WSI ♦»( I.S^TBlIAl-.^'^TBilA-n-TINF)/! 7^4 . ^♦'Wk^^DR ) 
WRITE(6»10)AI( IA),Bt( IA),CI( IA),DI( lA) 

CALCULATION OF TA 


I A=1 

BETA( lAj^Bl ( lA) 

CAMMA( I A)=Ol ( 1 A)/B1 ( I A) 

00 380 I A=2, 76, 1 

BFTA( I A) = B! ( I A)-( Al ( I A)*»CI ( lA-1 ))/BETA( ! A- 1 ) 
GANMA( I A) = ( Dl ( I A)-AI ( IA)*GAMMA{ IA-1 ) )/BETA( I A) 
330 CONI INUE 

TA{ /6)=GAMMA( 76) 

00 0 30 . 75, 1 

I A= 76-M 

TA( IA)^GAMMA( IA)-{CI ( IA)»TA( lA^I ) )/BETA( I A) 

U30 CONTINUE 


TEMPERATURE LIMI TER 


DO UUO 1 = 1 , 76. 1 
I F( TA( M.LT, 1000) 00 TO 440 
I I -* 1 

1 A( ! )- 1000 
440 CONTINUE 


CALCULATION OF TOTAL CURRENT 



I A= 1 


I =0 

I TOT=0 
TJ-TA( lA) 


Figure A-3. Continued. 
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CALCULATION OF CURRENTUENS I TY 

CALL CUROEN 
lT0T=(CD*^PI«^0R**2)/4 

00 450 1 = 1,75, 1 

1 A= I + 1 
TJ=TA( lA) 

CALCULATION OF CURRENTOENS I TY 
CALL CUROEN 

I TOT= I TOT'»*CD* I *P I ♦^2*DR*^2 
450 CONTINUE 

TC(L) = TA( 1 ) 

TT( L)=IT0T 
TV( L)=V 

CURRENT LIM'TER 


0=1 .000U0001*IMAX 

E=. 99999999* I MAX 

I F ( I TOT.GT.O) GO TO 490 

IF ( I TOT. LT. E. AND. I F. EQ. 1 ) GO 

GO TO 500 

490 V=( V*IMAX)/ITOT 
IF =1 

CO TO 500 

495 V=(V*IMAX)/ITOT 

If (V.LT.VMAX) CO TO 500 
V=VNAX 
I f =0 

500 CONT INUE 
IM= IM>1 

DO 540 1=1,76,1 
TB( I ) = TA( 1 ) 

540 CONTINUE 
RETURN 


Figure A-3. 


TO 495 


Continued. 
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Figure A-4 










APPENDIX B. FLOW CHARTS AND PROGRAMS FOR THE 
MICROCOMPUTER CONTROLLED, CONSTANT TEMPERATURE, 


SHORT INTERVAL IV MEASUREMENT SYSTEM 
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Figure B 



I. Get Cell Number Procedure 
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Figure B-2. Step Through Procedure 






Start Mode 0 


Figure B 
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4. Compute Current Procedure Mode 0 
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Start Mode 1 


Figure P 



-5. Compute Current Procedure Mode 1 




Figure B-6. Compute Current Procedure Mode 2 
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Figure B-7. Scope Procedure 
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Figure B-8. Print Procedure 
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Figure B-11. Tbout Procedure 
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/• flfilN fROGRflM flODUL' V 

/• 

nftlNFFFvGRflWROOULE. 

»/ 

r»£ fRGCFC-tjfiE FATSNflL 

cNO rvf; 

•j£T.:a: FSOCffiyRE fXTSM- 
cNO liTCHi 
ID 

ifOCFDiJRt E<TEm. 

END ID; 

:jcriCFU.."?»on. 

?POCfD<jfiE E<TERta. 

END 'jcTICai Ifm 
iTF?*MO(JGH; 

PFiXEDufiE E^NPL 
END STFPJTHROUGH; 

ECOPE- 

PFOCEDtJRF S<TE?.fa- 
END jCOPL 

PNFIX. PP-XE't^Pe E<TFfiM>L 
END P»,t 

POUT SROCcDtJRf EXTeSNflL 
ENC- *OUT; 

CECL^ :'€aS PVTE EsTEFNtiL 
DrCi fiRE C»m SVTS £<T£E"PL 
■>'CW?E T 6VTL 
EEuiN . CSu. jET'Ifli 
M CSSf COHffiND. 

-jO to ICfNT. 
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« 'MNTF.;-. ■<*:•! " =::: 

X’^vT'^RTm.' • :tfN 
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•:«.;£CT •Kxt.' : rf fCFl' IH F! CW, 
CCrlflLfi >OK£y ?«■ riMSS rj .Cflt*. iRC 




-3ET ixnm ftXnif 


/* 


COWfiW«*«CCHf; 


DO; 



i 1 «C4!W f.vTg R<TEFM1. 

: 1 f»fC4!»ft£ ?VTE 

4 1 SPIHT pC0CFDliK<UXfiTl0fi, WinCim 

Z 1 (.€Ci!«£ t.OCATICIN HftiFFSh 

i 2 cf oflRE fimm ?m 

7 2 £» WlHTi 

$ J ffTCH: PIWCFC4J(«*|j3CfinOH,;4UKh«5> rATtim^ 

? 2 DfCuVe LOaTIOM 4W*fSS» 

le 2 D€C1 liSf mrom SVTc/ 

n 2 ao FFTCJi 

'J 1 ‘ST01: ?*•« IC. 

3- 2 r-a«£ Cx*> 5VTE 

14 2 START CflLI. W;f»T<.C«’d>,3a>. 

15 2 cflu «tch 4 . r*. 2‘j 

16 ? cai ?rint\.':slp, 2'-. 

1* 2 IF co?ww><3> ■; leH nei 

IS 2 '30 TO STARTi 

Bi? 

’? 2 IF C0m»N&»3^ > :iH THEN 

26 2 130 TO START. 

BiF 

.'3 2 C0W«D'6> ' av^'6 - - !6n, 

22 2 fW 

22 1 ao -a)WlW<.*«C«.ir. 


*X^JL- Itf-IRWriON: 

■2{f. AREA •::? • «:•« tk* 

•'UlrfiF -REii ;I3 * ^C* 

:*Tii*:>. - '■^2:^ lo 

.'4 I !>ifr iff'' 

j ••Otfiirt r*RCF' : ■ 
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iS-tl ^ *1 . 1 COf^IwnTIOH *3f ?!OWl£ C£Ll IDftiC'U/ 

;KT ?«JiC© ift 

P**Lff IhVlj(v£& 8V' fL?l6<8 :F1.I&. £SC 
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11 I 

a l 
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14 1 
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CELLIIDmODUi. 

DO. 

oeClWE 97 SESSi2> byte FXTFFML 
D€OjiK LOTO eVTe eXTEWL 
DKliWE TVPf<2> 8VTF &TEFH(t; 

OeafifiE Ca.F<2) BVTE £<T«YflL 
WINT PWCEDOREtLOCATl'JfcMUICHRRS) EXTEPWAU 
DFCUWE lOOmON AODfiEaS; 

DECUiFE HUfVCHMS pVTE/ 

END rfilNT, 

FFTCh: ?«KE5XJW(LOC«Tm MUW>f»S) eXTSFJA; 
mM, uXATION 4DDFESS> 

DEOAFE NUN»CHMS BVTt 
EW rFTCHj 

/• 

/• 

ID 0POCEDUFE ?*JMC 

OECUWE !€«K«> BYTE MTA( «lOTl<4CELik' 
v« # ISC * W • » 

IH 4 TENP 

iiOJH Ti*' BVTE DATA. iIEli T>FF w » 
DfOj»£ t.*> BVTF WTft« LOT 4 ‘ • i; 

DfCL-tfE 9\*» BVTE WTm sTFESS ^ •; 

/♦ 

■3i TVPE 

, 41 

CfL, PFINT- T<-B>,^C-i». 
call ?ETCH\.TVfE.B),i>, 

4«IN 'XL 

IP •.rTfE*;*,' o ^ 

T-fE'O; » 

■jO 'v -tflN. 

iW- 

mm 

r.;c. 


V 


V 
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• • • • 
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^C.- .j’ ; 1 

!* i,jT ;• * ’iM' "’zi* 
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*• CWFitR 


i 
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41 2 

•>2 2 
■42 2 

44 2 

2 

■Hi Z 
*7 2 
3 

*} 2 

;« ; 


51 2 
22 2 

52 2 

54 2 

55 I 


L0T<1/ » lOT*2>' 

« TO iftST. 

EHD) 

lOT.2> ‘ «2«M. 

CAi. PRINT* CFLf-2>i 

. * 

^ * STFES3 liva 

r* 

CfiLl PRlNT<.S« 6 >»^rt>i 
CHJl FETOK. STilEsSdi)^ 1); 
THRO, cm PrTOH. STRESS* !).!)> 
IF STRESS*'!) O 20H TH& 

oa 

STRESS*'^/ • STSESSiDi 

Od TO TnRUi 

9i&> 

STRF3S\1> • 93m 
CaifRTNT<CRLF.2)i 

/• 

. « TRBi )€S&IN6 

/* 

Cm.MINT(.CRLe.2)i 
Cm.RR!N7i.nFr)0.7lD>^ 
cm PRINTS, a&f, 2); 

E» !0i 

END CaillOlTlMf. 
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rt_. CuilrlLr? 


ii 


IS-II vI 1 CCfiPILriTION \f ilODlii iTffTHftOlJGHfiCL'OiJ 

JsCT Pt«I5‘ IN Fi:57EP OBJ 

iNVOK® 3V: FlJlSB iFliSTEs SSC 


i ! 


' i 
. 1 


/♦ 


SiF THROUGH Htt-ULE 


♦/ 


: STEf^rHSOUGWHOim 

L*0/ 

J 1 DEQFlfiE ?0INTu55) STRUCTURE^ 

CWfRENT«’2SD.) IMPRESS) EI-ITERNflU 
: 1 DECLRfE SUM 3D0RE35 EXTERNMj 

4 1 DECLARE INCREMENT BVTE EXTERNAL 

5 1 OeCUiRE <LT,«0E.H) BVTE EXTERNAL 

6 1 D£CJ«E N BVTE EXTERNAL' 

r I DECLARE MEAS 5VTE EXTERNAL 

3 1 DECLARE NSAVE(25D) BVTE EXTERNAL 

i 1 DEaARE fCASUREHBff'lSD) STRUCTiJBE',’ 

P0t«S\4) BVTE) FXTERNAL 

•/ 

IB I OELAV PROCEDURE EXTERNAL 

U i 94D DQ.AV> 

.•♦ ♦/ 

L i PMAX; PfiOCET-URE EXTERNAL 

L* i Ef© PflAXi 

/* V 

14 1 TBOUT. 

pROCOURE ''iflLUE) EXTERNAL 

15 2 DEClPftE ^UE ADDRESS; 

iP 2 END TfiOUT, 

/* *t 

ir 1 CGhPUTEICURRENT: 

PROCEDURE E^TERNAL 
12 2 END 'lOMPUTcfaRREHL 

,• * 

iS 1 STSfTNRijWiH 

rfiOCEC'i.SE 

DEC* ARE A‘irAQ£*fOINTE* AC’i’RESSj 
DECLASS .ITC-E EASED -CtLTAC-SfPO INTER AODSL'S, 
•-ECLASS •TAPT.»'ICNV PVTS AT'0f""h.‘: 

DECLARE -TATijSIPOINTF? AC'CfESS. 

DECLARE STSTijf 8(5f£,*i :?AT'j£.*fOlNTE? -tTE, 
DECIARF JOS EVrL 
DECLARE ;T'.'P -vrE AT'C'E^Ah;; 

.’EC! ARE iNSEl PVTs AT'sFr^fn,'; 
lECLAPE 'V'SCR SV^S AT'tfF"r"A'. 

CELL ARE ICiijNTER RvTE. 

IECIApE "i-vREi; 

It * • 

‘-Erf* 


186 


.- 1 . 


12 


Cfil T50UT^PRtVI.3US>i 


^0 

^1 

:s 

•*6 

41 

* 2 . 

4S 

•W 

45 

46 
• 1 ' 
48 
4? 
38 

31 

32 

33 

54 


Moof = >jeHi 
IT : 86Hj 
ii s 88H; 
iurt * 68hi 
SfTiJP « m 
ijhScL '• 82Hi 

vOLTfWElPOINTER = eFTFDH; 

STflTbSiPOINTER 8FFFCH; 

INCREMENT = m 

N=8i 

POS « 8; 

DO rtHiLf POS » 8i 
CPU. CCffUTEfCURRENTi 
iXL DflPV; 
nUXPCR • dOHi 

R&iO: STflRTiCONV » 8FFHi 

M8IT: 

DO WHILf'STflTUS PNO i6M* ^ 9i 
EM) NAIL 

POINT(N). vOLTf«E(f€ft£) •- ilTGE W 87FFH; 

POS ^ SHR^VLTGE 8ND 80OH,11D1) 

IF ^(PREVIOUS PND 0FFEH)‘>?OINT<M>. yOl.Tf«E*;f€PS) W •.frrri" THEN 


Z'sf 

56 

4 


yui 

PREVIOUS •- POINT(N> ^X-TPiiE' 

57 

4 


IF COUNTER > 18 H I)€N 

Zd 

•4 


(jO to QUIT; 




FL«F 

55 

4 


COUNTER :■ COUNTER ♦ 1; 

a 



tjO TO REnD; 

a 



END; 

6Z 

“ 

i4jn 

counter ' 88H; 




IF node O 82 TNEN 

i< 



CPU. ?:«o 

65 



:4 ^ N -r i; 

66 



PREVIOUS = eeH; 

67 


END. 
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?<£nVtvii£rS/s^ 




■Xl TgOtjTvPREVIOUS); 
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:*TSr#Ti«ROUGii 



LNC' 

:7E?«TNSa''jh^N0DUli. 


r : 

ftfOSf'STIuN 
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r ;II 

's = r7'*D 


«(«« 

I-tL- -^E: 

^ fir? * 



6Tr*I: 

> fllr = -C’ 



• ‘ ^Hr{ 
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'*} - 

£?r 

' V : ‘ 
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ii-ii -I J i:0«FfL'»TIi?N OF HOf-IJli jcTCflLyWCWti 

JsCT fHiixtz .^!iC£r‘ Iff ?1. 3C£i.;_ 'i6J 
I'l'irlLTR IffVOf.ti* ?V. ?tf)86 fi.'jCPu.. SRC 


ijcTlCFLLMJfPflOr'Ux. 


PRINT rROCE-UFEvLOCanOH. M.•HC«iRS:•E^ra^ff^iL 
OECLiiSE LOCfiTION POOftfSS, 

DEaPRE mJHCitiRi 3VTE.. 

END PRINT. 


rETCH .PR'OC£ri«E«xuCRTION. HUftfCrtifi: .'i‘.T£RNHL 
if'XhRE tOCliTION flC-OfE-5. 

DEa«E NUHFCHflfiS RVTF, 

•fjO i^CN; 


CECiflftE CPIF«.2) SVTE S<T£:*«L 
CfOliRE CELLOS) PVTE E<T£5NftL 
uECLRRE I PVTL 

C€CLffi£ At*) PVTE MTAvIELi W-lfftFR - s 
DECUWE Tv*) PVTE DiiTA v'tm TEHP » .« 

DEaiRE TEflPvG) PVTE EtTERNfiL 
ijET*CELlfNW1- 

PROC££-URE PUBlICi 

CFiX'.0.‘.Cai(i).CEUv2).Caj-vS).fFLl',4> = 020H; 
I = 0; 

CfiU PRINTt :v0>..0EH); 

'IRLl rFTCHv. CELLvd*-*); 

IF 'CElLvi.' * 0C4i» THEN 
Wi 

lEil - CftL'0)' 

CE;,L'.0).CELl'i) ■■■ 020H. 

'iO TO PET. 

END; 

'Iru PETCH* . CELi.' 2 ) • J i 
IF '.CELl'.I) - OC-H.' THEN 
I"); 

'IEIl' 2'' * CFuvl.'. 

■IEj-vI; = 'IS.:.' 0 •. 

CEii.'.P' • il'jhi 
•kT *0 JET. 

-IF. ‘I-ll ~::-v 
IF . II 'I'Of’,* 
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->3 i cM:* ‘£iKFu fM.'M?WC‘ULr. 
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? fifISLr nKrn SIZF = ■?Wlrt 
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3e LitfS FfnD 
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il i 'A-«PIiilTION -jF flODiJlF CiWjTSCUPS&TK'WJ 
.ipJECT >!0Cn? *•. Ih Fl.CCXfiT OBJ 
^if?Iixr INVOfEi* ?.y; PL.189 .FI CCRHT. ESC 


" «/ 

*■/ 


C4JHPUTE OiPFENT HOCtlLi 


/% 

,-k 

, MwvwiWiitiiliMntrM..: • 

J CGf!P'j 7 E«Cl^RENTIWDtA£. 

i)0; 

. * 

: J 'piEoT 6 F 0 CF<X«E (VfiuE) E^ 7 t?f(fiL 

; OtaiiSE yft. uF SCi-RESS/ 

<1 ;■ thO TK*JT.i 

K 

3 i DfCUfiE F 0 IhT' 255 > STRUCTURE- 

y 0 lTaG£- 25 Dy rtCi-REsS, 
ajR 8 £? 4 T< 2 ;ii. tiCC-RESSJ E>!TERNftL 
0 1 CfCiJiRE flE^S 6 VTE EXTEWiLj 

r i CCCLiWF vLLnCDtH) FVTE EsTERNfiL 

0 J DEasiRE INCRSfNT BVTE EXTERNfiL 

1 i OECliiRE sun liDDRESS E-iTERNPL; 

li i DEC.PRE H SVTE EXTERNAL 

ii 1 ifiiiV FROCEt'UR? rUBLIC, 

2 DO; 

;l‘ : DFCLflRE P BVTEi 

i-i : M aHILE *.P< 0 FFH)i 

i 3 •) P = ? ♦ L 

-9 4 EN 0 < 

i 7 ■ ENOi 

IS 3i> i-fliiV; 

L- 1 :GHP‘.iTEs:UPRENT: 

SROCET-ORF Pi-aiC, 

ii 1 KCbiRE ',niN.Mft>i> BVTE, 

;i ;■ DEilPRE fiWE *CR£S 5 > 

/ ■'ECL'iRE CELTS RVTE; 


*/ 


DO csss i-lCt'E, 

.'ftycO M; 

Ir Ll ■■ i rr€N 


rOi*l 7 • ,N,' Cl'PPENT ' !•!£•£ - 

iv.1 ;v:i • 15H. 

ENi'. 

E:..-E 

*jvt - •’Ih; 

•vii ■• run -i*>. 

■y, 

i-j ■'■'XFi 

Em* ‘ 0 ; -i 3 . 
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•1. ^50 CuflflLiR 


li 


iT0C£i DO; 


,1 


I? lT : 0 rh£N 

41 

4 


DO; 


0 


iwi s am -oiH, 


s 


POINTvNa CUPfiaT-flOiS/ s SW!; 

44 

s 


c«ll rBooT\ai«j, 

4C 

s 


END; 




ELSE 

4b 

4 


DO; 


3 


NODE • 02H; 

,40 

%• 


SDH: SflVL- 

49 

3 


60 TO nOOf E; 

v0 

5 


m 

5i 

4 


END flODfl; 

ci 


itJCf;; DO; 

SI 

4 

HIM 

» 62Hj 

54 

4 

nix 

- m 

55 

4 


IF rt<i 1 n€N 

:i5 

4 


DO CftiE Hi 

57 

V 

/ ^ 

DO; 

50 

0 


slifl ‘ SWI ♦ 01; 

59 

0 


FOINT(N), CURRENKHEre) = SWi 

a 

0 


CftU. TB0UT<5UM); 

a 

0 


END; 

0< 

J 


DO; 

il 

0 


SOM am 04H; 

04 

0 


POINT* N),CURREHT^f€ftS) = SOft 

03 

0 


Cfll TK*UTt5U1)i 

00 



END; 

07 

‘J 


END; ;• ENDCflSE */ 
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IF lOaTii > NIN) PND DaTfi •; m) Ti€N 

n 

5 


Da 


0 


SUN * aJH ♦ INCRE«IT, 

75 
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EtDlNT^Ni iaSRENT^REPi) - fwT; 

*4 
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CflLL *kWai«;; 
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END; 
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IF 'rtlTfl -‘IN' *r€i 
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-<> CvrlfLr? 


iiihll c.: :f “iX-.itf KijriKw^ 

■:.r;--:* .‘iODwf it-iCfy in FJ rfCCri.jS; 

CCmFIlF? iNrW,Et* tv fi.iis« FJ.?f::<F£ ir: 


iiFCUiofCOPF £*ia?t!ir 


KOPfWGyOLF 

OO; 


: i Cl FFOCS-fJRS 3VTF £NTr?Nnu 

3 <’ FNO Cl; 

4 1 PLOT: FROC^cRf'.nXIS'* rI\T?Sf^; 

5 2 OfCLNfiE *IS fVTF; 

i ? FND PLOT; 

/* 

7 1 DECLflftc rOlNTuSG'* STPUCTi.«Et 

V(XTfiG£<.2:3>,. (W<?f55; 

CUPR£NT(2!»> ii0f-Rf?5> F^TFnA; 

9 1 DfCLflSE NSi9VEv2Si 8VTE EXTERNAL. 

9 1 CfCLnRE MERE ttVTc EXTESWb 

10 1 SCOPE. 

PPOCFT-tjPF P(.«l fC; 

11 2 Cfa«E I eVTE, 

12 2 t)Ea*E PIT PVTE; 

IS ; OfCbiPE »OflC2La.C-)iCiLO> Fff-SFSs nT*0Fr?.>h..': 

14 2 OO uINPuT\(j€Cn> liND 02h j ■ 

15 3 itfE?' 

M I -0 TO •NSnVE'?€0S’-2-: 
le 4 IiiCllO ■- POINM,'. 'OLTP*x-fSS-i 

ir 4 C-PC2L0 " POIN^^^afPfNT•;■ff•r^■. 

18 4 rM> FhEFF; 

1? I ENfii 

20 2 PIT ‘ Cl 

21 ; IF rl T » ?■ r^N 

22 .■ CPU. r1 OTvNtPS>i 

22 2 tK> -CCff. 

24 1 :N() iCCFE^IHiUF; 




.^LT • ;^r 

-7K- 


r'^L’ jf ^ 


-•w fh*6 CC'Hf IlB? 


iiii-tl v3. 1 C4»!LiiTI0N OF flOCtle 'HTFSIWt€?S55WOOLH 

?lCv.n5 ?ljiC?y IN .rl.0WR.i56J 
'VjNriiJr I.Hv'OKeD 6V. FLN9J .Fi.ijPliR.sSC 
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J 1 


i i 
4 1 

: 1 
i 1 
7 } 

« i 
• 1 
16 1 
U 1 
H 1 
U 1 
14 1 

:5 j 
le 1 

1* J 
16 1 
1 ? 1 
1? 1 
a 1 
a 1 


1 


ij£TIPl«f^€TeRs‘mjDiJL£. 

uO; 

DFCLPRE P0INT(255> STFUCTUREv 
V0LTA0£<2S» RDORESS. 

CURR£NT^;^D) «ORESS) £<TFRNRL 
DECUWE LT SyTE E<TPfaib 
DECLPRE 'm(5> sVTE cXTtPNfiL 
C€CU«£ TPPEICH BVTF EXTKNiL 
OfCUWe TRP£K3FF 3VTE EXTERNAL 
C€aiiR£ STRESS^:) BYTE EXTERNAL 
DFOifiE LOT^S) SVTE EXTERNAL 
CECUWE TVPE<£) SVTF EXTERNAL 
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APPENDIX C. ANTIREFLECTIVE COATING AND 
COLOR MEASUREMENT THEORY 
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APPENDIX C. ANTIREFLECTIVE COATING AND COLOR MEASUREMENT THEORY 
Antlreflectlve Coating Theory 

Before sunlight reaches the energy producing area of a terrestrial solar 
cell. It mist pass through several optical media. Each boundary between media 
with different Indices of refraction causes the light to experience 
reflections which are non-uniform over all angles of Incidence and 
wavelengths. For unencapsulated cells, the media with the most significant 
effects are: 

Atmosphere 

Anit-reflectlve films on the cell's surface 
Solar cell material 

Sunlight which strikes the cell's front surface Is partly reflected and 
partly absorbed as shown In Figure C.l. An antl-reflectlve film is a thin 
coating, applied by some manufacturers, to enable the cell to capture a 
larger percentage of the available energy. Ideally, the antl-reflective (AR) 
coating should reduce the reflection loss to a maximum of three percent, and 
result In a greater electrical output by the cell. 

A typical AR coating Is optimized for the wave length of the sun's peak 

spectral response (C.l). A thin film of magnesium flourlde (M9F2) or 

similar material Is vaccum deposited on the surface. This film should have a 

thickness which Is an odd multiple of X wnere A Is the wavelength In the 

AR coating. The effectiveness of the coalng Is determined by the angle of 

Incidence at which sunlight strikes the cell, the Indexes of refraction of the 

atr and the film, and the amount of the spectrum over which reflection Is to 

be minimized. For light passing from air Into the AR coating at a right 

angle, the reflectance, p , Is given by: 

2 

7 
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Figure C.l Reflected and Transmitted Rays at the Surface of a Solar Cell. 
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where Is the Index of refraction of air and n2 Is the Index of 
refraction of the coating. The value of p is the portion of light that is 
lost in the transition. If the angle of Incidence is not normal to the cell 
surface, then p is modified by Fresnel's formula. 


2 1 

tan 


• 2 , 
sin 

K - ^2' 

2tar.^ 1 

k *2. 

2 J 

2s in 1 

h ^ * 2 , 


The angle the angle of incidence at the boundary and 4*2 is the angle of 

the refracted wave with the normal to the cell surface. These angles are 
related by Snell's Law: 


n^ sin())^ ■ n2 sin4i2 

The incident ray which passes through the coating is then partly reflected 
again at the silicon surface. This reflected wave then returns through the AR 
coating to combine with the original reflected wave (Figure C. 2 ). If the path 
difference between the distances that the two waves travelled is exactly nX/ 2 , 
where n is any odd integer, the waves cancel each other. 
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Figure C.2 Tnterfe *nce by Reflection from a Thin Film, 
Assuming an Extended Source. 


! 

Color Theory 

Images are seen by stimulating the retina of the eye with light energy 
in the 400 to 700 nm wavelength region. Listed below are the colors of the 
spectrum seen by the eye at various wavelengths. 


Violet 

400-23C 

nm 

Blue 

430-485 

nm 

Green 

485-560 

nm 

Yellow 

560-585 

nm 

Orange 

585-610 

nm 

Red 

610-700 

nm 


In order to categorize one object from another by color, certain 
differences must be observed. These are hue, lightness, and saturation (C.2). 
Hue is the exact shade of color seen and is a linear function of the dominant 
wavelength reflected by an object. Lightness deals with lighter and darker 
shades of the same hue. It can range from a perfect white to black. 

Saturation pertains to thf. pureness of a color. It measures the degree of 
difference of a color from gray of the same ]*ghtness. The interrelationship 
of these three quantities can be visualized by using the diagram of Figure 
C.3. 

Also ImpOvtart in color perception is the spectral distribution of the 
light source. The same object may appear a different color when viewed under 
a different type of illumination. The spectral characteristics of two common 
illuminants are shown in Figure C.4. Illuralnant A approximates the spectral 
distribution produced by a tungsten lamp (indoor illumination). Ii- produces 
greater energy at longer wavelengths, particularly In the infrared range. 
Illuminant C, however, has a more nearly flat response with a small peak at 
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about 450 nm. This more nearly approximates the solar spectrum (outdoor 
illumination) in the visible region and produces a slightly bluish hue. The 
total solar spectrum is shown in Figure C.5 and it can be seen that only a 
limited portion is useful for photovoltaic conversion, and only a limited 
portion of that, can be seen with the eye. 

Measurement of perceived color may be expressed in severa? different 
ways. The most straightforward approach is to state the percen' spectral 
reflectance. This may be presented either in tabular form listing percent 
reflectance for specific wavelengths, or as a continuous plot over a 
particular range of wavelengths. 

Another more advanced form allows the apparent color to be expressed as 
a combination of the three colors of light. These are known as tristimulus 
values: x for red, y for green, and z for blue. Each of these parameters 
range from 0 to 100; the higher the value, the more intense that particular 
component is. A similar representation is by chromasticity values. These 
coordinates, x(red) and y(green), give the per unit fractions of these two 
colors reflected by the smple. However, no information on light intensity is 
obtained . 

Another widely used representation involves dor space coordinates. The 
roost recent method accepted by the International Commission on Illumination is 
the 1976 CIE L*A*B formulas. The L*A*B values tnat are obtained can be 
plotted in three dimensions, as shown in Figure C.6, which aids in determining 
color differences. The L parameter is the overall lightness of the sample and 
has a range of 0 to 100 relative to the measuring device’s calibration 
standard. The A value reveals the degree of green or red reflected by the 
sample. A green hue would result in a large negative A, while red is a large 




Figure C.5 Spectral Distribution of vSolar Radiation 


L 

white 



Figure C.6 197b CIE L*A*B Color Space 
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positive A. Similarly, B is a meausre of the blueness or yellowness 
reflected, with blue in the negative scale and yellow on the positive scale. 

Eich type of color representaton can be obtained from the IBM measuring 
system. The 7409 unit divides the visible portion of the spectrum into 31 
channels each of 10 nanometers width. The fiber optic detector divides the 
illuminated region of the sample spatially. During measurement each spatial 
region is samples over the 31 frequency regions using a 12-bit A/D converter. 
The sum of the different spatial readings for each frequency interval is 
stored. After all measurements are taken, this sum is divided by the number 
of spatial points measured to provide an average value for each frequency 
channel. The voltage V(I,J), for wavelength I at measurement position J, is 
linearly related to the porcent reflectance R(I,J). The constants in this 
relation are determined during the device calibration procedure. 

In order to calculate tristimulus values x, y, and x, the Steam’s 
Coefficients T (I), T (I), and T (I) must first be determined. They are found 
for each wavelength by multiplying the spectral trlstiimilus values of an equal 
energy spectrum for a particular observer angle by the relative spectral 
irradiance, E, for a specific illuminant. Tables, listing these numbers, are 
published by the National Bureau of Standards. The following example computes 
the X Steam’s coefficient at 500 nm for illuminant A with a 10^ observer 
angle. 

T, (500) = (500) X E, (500) 

X lU A 

- 0.0038 X 59. 8t. = 0.22747 

The percent reflectances and the Steam's Coefficients are then used to 
determine the tristinrulus values using the formulas given: 
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I 

i x(J) » (T (I) X R(I,J) X Cl/ T (D) 

I X j X 

y(J) » (T (I) X R(I,J) X C2/ T (1)) 

I y I y 

z(J) « (T (I) X R(I,J) X C3/ T,(D) 

I ^ I ^ 

Cl - 9G, C2 • 100, ami C3 - 116 

Chromas deity values, Cx(J) and Cy(J) may be realized from the relations: 

Cx(J) - x(J)/(xU) + y(J) + z(J) 

Cy(J) - y(J)/(x(J) + y(J) + z(J) 

The 1976 CIE L*A*B color space coordinates are also found from the tristimulus 
values . 

L(J) - 116 X (Cx(J)/10000) ^^^-16 

A(J) = 500 X (Cx(J)/9800)^^^-(Cy(J)/1000)^^^ 

B(J) = 220 X (Cy(J)/10000)^'^^-(Cz(J)/11600)^^^ 

These parameters are used in describing the observed color changes of solar 

cells- 

References 

C-1 Solar Cell Array Design Handbook p. 246-248. 
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APPENDIX D. COMPUTER PROGRAM LISTINGS FOR MODELLING 


ENCAPSULATED AND UNENCAPSULATED CELLS 
IN THE SECOND QUADRANT 


PRECEDING PAGE BUVNK NOT FUMED 



//RAHRTB JOB (0915-2-992-00- , 18 : 00, BOO ) / ROB HARTMAN* 

// EXEC FTGICIG 
//C.SYSIN 00 * 

(;;*##*#*###<»#######♦« ####•»##*###*########## !♦###*###* ##*####<♦*# 

c** BARE ** 

C*^ THIS PROGRAM CALCULATES THE TEMPERATURE PROFILE OF A BACKWARD- 

BIASED SILICON SOLAR CELL, THAT IS NOT ENCAPSULATED. ** 

A GAUSSIAN DISTURBANCE WITH DIAMETER DEL AND AMPLITUDE LF IS A 
C*^» PARAMETER THAT CAN FIT THE MODEL TO A MEASURED l-V CURVE. ** 

Q»* ## 

Q»# «#*#«####«#«« 

C** MULTIPLE RUN CONTROL 

#####«##»##«»«###•»#««####»«»##««»####### 

c#**#*#*###*#*#*###*^***»*#####*»##^#***##*#*##*##*#*#^##***##*#**#*####* 

INTEGER I . IF.K.L.M. lA. IL, IM 
REAL TJ,WK,CO,DR,V, I TOT , TB, TA, VMAX, TC, TT 
C. IMAX.OT.X.AI .81 .Cl .01 . BETA, GAMMA. 0. E 
C,HL, ;RR. 1 1NF,WSI ,PI ,C,A.OEL,LF, TV, VAR 
COMMON I, IF.K.L.M. lA.O.E. IL.A. IM 
C.TJ.WK.CD.DR.V. ITOT,TB( 101).TA( 101 ), VMAX. TC{ 100).TT( 100) 

C. IMAX.OT.X.AI ( 100),BI( 100).CI( 100).DI( iOO),BETA( 100),GAMMA( 100) 
C.HL, IRR.TINF.WSI . PI ,C,DEL.LF.TV( 100). VAR 
DEL=25 
WSI=.04 
VAR=1 .2 
LF=0.0 
CALL MN 
LF=0.2 
CALL MN 
LF=0.4 
CALL MN 
LF=0.6 
CALL MN 
END 

SUBROUTINE MN 


##«###««#########»#### 4 #######»##############»#### ############## 

MAIN 

^^#########«# ###«################################««##*##«« ############### 

I NTEGER I , I r . K. L.M, I A. I L. IM 
Rf AL TJ.Wk.CO.DR. V. ITOT. FB. TA. VMAX.TC. TT 
C. IMAX.OT.X.AI .01 .Cl ,DI . BE I A, GAMMA. 0. E 
C.HL, IRR.TINF.WSI.PI . C . A. DEL . LF , T V. VAR 
COMMON I. IF.K.L.M, lA.O.E. IL. A, IM 
C. TJ.WK.CO.DR. V. I TOT. T8( 101 ) , TA( 101 ) , VMAX, TC( 100) . TT( 100) 

C. IMAX. DT.X. Al ( 100) .81 ( 100) . Cl ( 100) , Dl ( 100) . BETA( lOO) ,CAMMA( 100) 
C.HL. IRR.TINF.WSI .PI .C. DEL, LF, TV( 100), VAR 
IM=0 
HL--.003 
IRR=. 1 
TINF^25.0 
PI--3. 141592634 
A- 1 
I F^O 


I L = 0 

OIU .0S9602MI9 
IRR/HL+T INF 
no h 1=1.76.1 
TA( I )=C 
1B( I )=C 
5 CONMNUE 

DO 6 1=77. 101 . 1 
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TA( I ) = r INF 
TB( I )-T INF 
6 CONTINUE 

v=n.o 

VMAX=50.0 


CALCULATION OF TOTAL CURRENT 




IA=1 
1=0 

I T0T=0 
TJ = TA( 1 A) 

CALL CURDEN 
i T0T=(CD^PI*UR**2)/U 
DO 160 1=1,75. 1 
IA=1+1 
TJ = TA( I A) 

CALL CURDEN 

I TOT= I TOT+CD* 1 I *2^»0R^*2 

160 CONTINUE 
CALL OUTPUT 
V=25.0 
1MAX=1 .2 

RUNNING CONTROL 

C 

DO 200 K=1 ,40, 1 
I MAX= I MAX+ . 1 
CALL RUN 
CALL FINE 
CALL OUTPUT 
200 CONTINUE 
RETURN 
END 

SUBROUTINE RUN 

('I #-#««##################«####### H ###### 

C^^ NEW TEMPFRATURE RUN AND STEPSIZINC ** 

(^ffffff############ ###ff*lf ## If #ff 4## ##ff##ff#ff ff#ffff# ######## *#ff ##*#### ffffffffff 
Q If # # ff ff # ff M # # # # # ff ff ff # -If # ff If # ff » # # # # If 4f # # # # # # -4 # ff ff # ff ff If ff ff ff ff ff # # # ■» 4f # # # ff # -If ff # # ff ff # # # # -if 4f 

INTEGER I . I F . K. L.M. lA. IL. IM 
REAL TJ.WK.CD.OR.V, I TOT, TB, TA, VHAX, TC, TT 
C. IHAX.Dl ,K.AI ,BI ,CI ,01 ,BF TA,GANMA.D, 1 
C.Hl . IRR,TINF,WSI , PI ,C,A,nFL,LF. IV. VAR 
COMMON l.lF.K.L.M.lA.O.r.lL.A.lM 
C, TJ , WK.CD. DR, V. I TO F , T B( 1 0 1) , T A( 101 ) , VMAX, TC( H)0 ) . TT( 100 ) 

C. 1MAX.DT,X,AI( 100), Bl( 100),C1( 100), Dl( 100),BETA( 100).GAMMA( 100) 
C.HL. 1 RR, Tl NF.WSI , PI . C . DEL , LF . T V( 100 ) .VAR 
DO 161 L=1 ,20, 1 
D1 '2 

CALL NFWTEM 

161 CONllNUL 

DO 16P L=21 , 40. 1 
DT=1 

CALL NEWTEM 

162 CONTINUE 

DO 16 J L=41 ,60. 1 
DT- .« 

CALI NLWIEM 

163 CONTINUf 

DO 164 L=61,rtO,l 
DT - .6 

CALL NEWTEM 

164 CONTINUE 


DO 165 L=81 , 100. 1 

or-. 5 

CALL NFWTEM 

165 CONTINUE 

DO 166 L=1 ,20, 1 
DT=.4 

CALL NEWTEM 

166 CONTINUE 

DO 167 L=21,40,1 
DT=.3 

CALL NEWTEM 

167 CONTINUE 

DO 168 L=4 1,60,1 
DT=.2 

CALL NEWTEM 

168 CONTINUE 

DO 169 L=61 ,80, 1 
0T=. 1 

CALL NEWTEM 

169 CONTINUE 

DO 170 L=81 , 100. 1 

DT=.08 

CALL NEWTEM 

170 CONTINUE 

DO 175 1=1 . 100, 1 
TT( I ) = 0 
TC( I )=0 
TV( I 1=0 
175 CONTINUE 

DO 181 L=1,20,1 

DT=.06 

CALL NEWTEM 

181 CONTINUE 

DO 182 L=21 ,40. 1 

DT=.U4 

CALL NEWTEM 

182 CONTINUE 

DO 183 L=41 ,60, 1 

DT=.02 

CALL NEWTEM 

183 CONI INUE 

DO 184 L=6 1,80,1 

01^.01 

CAN NEWTEM 

184 CONt INUE 

00 185 i.= 81 , UK). 1 


DT= .01 
CAL'- NLWiEM 
185 CONTINUE 
Rl I URN 
END 

SUBROUTINE ^INE 

C**n*»»*******» '••***'*'•****•**'•***** ****" ****<******»***'►'***'•'•'*******'**•**'•'**•** 
(;*#*********###•»#**■»•»#****»•*'*'< *1 '♦*'♦** *♦'♦*********'***** ******************** 

EMRA riNr STEPPING ROUTINE F OK BETTER PRtCiSION 


INTEGER I. ir.K.l ,M. lA. U . IM 
REAL TJ.WK.CD.OR.V, I TOT. TB, TA.VMA.X. TC. TT 
r, IMAX.DT.H.AI .TTI .0 1 .01 ,IM T A , GAMMA . 0 . E 
G, HE. IRR. T I NF . WS I , PI . G. A. on , LF . TV. VAR 
COMMON I , I F . K , I, . M . 1 A . 0 , L . I L , A . I M 
C. TJ.Wk.CD.OR.V. I TOT. TB( 101 ) . TA( 101 ) .VMAX.TC( 100) . TT( 100) 

C. I MAX. OT . X. A I { MU) ) , B I ( I 00 ) , C I { 100 ) , 01 ( 100 ) , BET A( 100 ) . GAMMA{ 100 ) 
C.HL. IRR. TINF.WSI . PI .C.OEL. I F. TV( 100) .VAR 
DO Ub 1 = 1 . UH). 1 




P(H*P ' AS 

^ A*pAi./'n 
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TTi 1 )=() 

TC( I )=() 

TV( I )=0 
175 CONTINUE 

DO 181 L=1 ,20, 1 
DT=.0n8 
CALL NEWTEM 
131 CONTINUE 

DO 182 L=21,40,l 

DI=.006 

CALL NEWTEM 

182 CONTINUE 

DO 183 L=41,60,1 

DT=.005 

CALL NEWTEM 

183 CONTINUE 

DO 184 L=6 1,80,1 

DT^.n03 

CALL NLWTEM 

184 CONTINUE 

DO 185 L=81, 100, 1 

DT=.001 

CALL NEWTEM 

185 CONTINUE 
RETURN 
END 

SUBROUTINE NEWTEM 

(;#**************•»»**#*•********»**«•»*********■*»*********************•***** 
c** NEW TEMPERATURE CALCULATION ** 

Q***»******« *****#*****(*»*****(*******************#******#*»*****»******* 

(-*»(MH********#<********»*»***#**»****#**********#***#***#******<‘«********* 

INTEGER I , IF.K.L.M. IA. IL, IM 
REAL rj.WK.CO.DR.V, I TOT , TB, TA. VMAX, TC, TT 
C. IMAX.OT,X,AI .81 ,CI ,DI . BETA. GAMMA. D. E 
C. ML. IRR.TINE.WSI.pl , C . A, DEL . LE . TV. VAR 
COMMON I.IE.K.L.M.IA.O.E.IL.A.IM 
C. TJ.WK.CD.DR.V. I TOT . TD( 101 ) . TA( 101 ) . VMAX. TC( 100) , TT( 100) 

C. IMAX.DT. X.AI ( 100) .Bl ( 100),CI ( 100) ,01 ( 100) ,BETA( 100) ,GAMMA( 100) 
C.Hl , IRR, T INE.WSI . PI .C.DEL.LF, rv( 100) , VAR 
UWORMATI’ *,'A B ’.FU).4.’ C ',F10.4.* D '.F10.4) 

I A=1 
1=0 

TJ=TB( IA) 

C 

C CALCULATION OF CURRFN T DLNS I T Y 

C 

CALL CUROLN 
C 

C CALCULAflON OF THERMAL CONDUCTION 

C 

WK=WSl/( .08+TJ/350) 

C 

C CALCULATION OF FHFr^MAL DIFFUSION 

C 

A=WK/( 1 .647*WSI ) 

C 

C CALCULATION Of FIRST SFT OF MATRIX FLEMENTS 

C 

AI( IA)r() 

Bl ( 1 A)^-2/0R'**»2-1/( A^DI ) 

{ I ( I A)^;VDR^*2 

01 I IA) = (2/DR^^P-1/( A^nr ) )**TJ-2**TB( I A-M )/0R*^2 
C- ( cn»^V-HL^T J+HL T 1 N M- 1 RR ) /W»v 
C WRI TE(6, 10)AI ( IA) ,BI ( IA) .Cl ( IA) .01 ( IA) 

DO 300 1-1,74.1 


:2b 


onoo onoooooooooo o ooonocoooooo 


IA=I+1 
TJ=TB( lA) 

CALCULATION OF CURRLNTDENS I TY 
CALL CURDEN 

CALCULATION OF THERMAL CONDUCTION 
WK=WSI/( .58+TJ/350) 

CALCULATION OF THERMAL DIFFUSION 
A>WK/( 1 .647*WSI ) 

CALCULATION OF THE REMAINDER OF THE MATRIX ELEMENTS 

Al I IA)=( 1-1 )/( (?*l-l )^DR^*?) 

Bl( IA) = -l/DR*^2-1/( A^DT ) 

CM IA)=l/( (2*1-1 )*DR**2) 

DM IA) = -AI ( IA)*TB( IA-1 )+( 1 /DR**2- 1 / { A*DT ) ) *TJ 
C-CI ( IA)*TB( IA+1 )-(CD*V-Hl*TJ+HL*TlNF+lRR)/WK 
C+( TB( IA+1 )-TB( IA-1 ) )**2/( (205+TJ )*4*0R**2) 
WRITE(6»10)AMIA),BI(IA),CI(IA),D!(IA) 

300 CONTINUE 
1 = 75 
IA=I+1 
TJ^TB( lA) 

CALCULATION OF CURRENTDENS I TY 
CALL CURDEN 

CALCULATION OF THERMAL CONDUCTION 
WK=WSI/( .58+TJ/350) 

CALCULATION OF THERMAL DIFFUSION 
A=Wk/( 1 ,647*WSI ) 

CAICUIATION OF TMF 76TH MATRIX El FMENTS 

AM I A) = 1/(2*0R**2) 

8I( IA)=-1/(2*DR**2)-1/( A^DT) 

CM I A)^0.0 

DM IA)^-AI( 1A)*1B( IA-1 )-(BI( )A)+?/( A*DT) )*TJ 

C-( CD*V-HL*rj*^Hl Nl ♦ I RR )/WK 

C+75. 5*Hl *WSI *( 1.5*^rn( IA)-.5*TB( I A- 1 ) - T I NF ) / ( 74.5*WK*DR) 
WRITE|6.10)AI( IA).BI( IA).CI( IA),DI(JA) 

CALCULATION OF TA 

IA-^1 

BLIA( I A)-BI( I A) 

CAMMA( IA)=DM IA)/BM I A) 

DO 380 I A-:2. 76. 1 

BFTA( I A)^BI( I A)-( Al( I A ) “Cl ( . A-1 ) )/PFrA( I A- 1 ) 

GAMMA( I A) = ( DM I A)-AI( I A)*GAMMA( I A-1 ) )/BETA( lA ) 

380 CONT I NUr 

TA( 76)=GAMMA( 76) 

DO 430 M-1 , 75, 1 
I 76-M 

TAI !A) -GAMMA( I A)-{C1 ( IA)*TA( I A+1 ) )/BriA( I A) 

430 CONTINUE 


an n nnn ooo oooo oo 


V 


TEMPERATURE LiMl TER 

DO 440 1=1,76,1 

IF( TA( I ).LT. 1UU0) GO TO 440 

I L = 1 

TA( I )=1000 
440 CONTINUE 


CALCULATION OF TOTAL CURRENT 




IA=1 

1=0 

ITOT=0 
TJ=TA( lA) 

CALCULATION OF CURRENTDENS I TY 

CALL CURDEN 
lTOT=(CD^PI^OR^*2)/4 
DO 450 1=1 , 75, 1 
I A= I +1 
TJ=TA( lA) 

calculation or CURRENTDENSITY 
CALL CURDEN 

I TOT^ I TOT+CD^I ♦^PI^2^0R^*2 
450 CONTINUE 

TC(L)=TA(1) 

TT(L)=ITOT 
TV( L)=V 

CURRENT LIMITER 

############ If ########«########### ######«############## 


0- 1 . oonoooo 1 * I MAX 

E=. 99099999* I MAX 

If ( I 10! .GT . D) GO fO 490 

IF ( I rOT . LT, E. AND. I F. EQ. 1 ) GO TO 495 

GO 10 SOO 

490 V(V*IMAX)/tTOT 
I F =1 

GO TO SOO 

495 V--( V*IMAX)/ I TOT 

I ! ( V. L! . VMAX) GO 10 500 
V=VMAX 
I f -o 

500 CONTINUE 
IM= IM+1 

DO 540 I =1 , 76, 1 
TD( I )=fA( I ) 

540 CON! INUE 
Rf TURN 
t NO 

^'UBROUTIN^ OUTPOr 

OUrPUl IN lAREL ! ORM 

(^lfl^##lf#lflflf#-if##*f«f#4»lflflflf«f#«f4f<<fM'lftf4f«ffflflf#^4lflf%4f»»lf###lflf>f«4f#lf##IMf#lflflflf#lflflf4f#ff# 

i N! I n R I . I I . I' , I , M, I A . I L , t M 
Rf AL rj.WK.Cn,OK. V. I lOI. TR. f A. VMAX. U\ M 
C, I MAX. 1)1 .X.AI .Bl .Cl .[)! ,B! I A , GAMMA , 0 . L 
C, Hl , I RR. T I NF . WSI . P I , C. A. 0! I , L F . ! V, VAR 
COMMON j.If.K.L.M, IA.0.L.IL,A,IM 
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C, TJ.WK.CD.OR.V. Mor, IU( 101 ), TAi • ' i ,VMAX. )C( 100 ), TT( 100) 

C, IMAX, OT.X, At ( 100) ,Bt ( 1C ‘) ,1, ’ . 01 ( 100 ) , BETA( 10O),GAMMA( 100) 

C,HL, IRR,nNF,WSI.PI,C,OEL,l' iiO,.VAR 

I F( IM. EQ.O) GO 10 1 
CO TO 9 

I WRI TE(6,2) IM.VAR 

C rORMAT( ' 1 ' ,I1X, ' TIME = . . • ' INCREMENTS, THE VAR I AT I ON' , 

O' I-R';M IHE standard (,;■ ■ »J ;; ',F12.4) 

WR. I ?.{'), 3 )DEL, U 

3 F 0 H? 1 A:i* '.Rv.’fHE C,f\ : I IS ' 

C.F12.4,' WIDE. ANO H,’. O AHFlITUOE OF ’,F12.4) 

CO TO 15 
9 WRITE(6, 10) IM 

10 FORMAi ( ' 1 ' ,4X, ' TIME - ',I12,2X, ' INCREMENTS') 

15 WRITE(6,20)V 

20 FORMAT( ' ' , 4X, ' VOL fA,t:C = ' , F12. 4, ' V. ' ) 

WRITE(6,30)ITOT 

30 FORMAT( ' ', 4X, ' CURRENT = ' , F 12 . 4 , ' A. ' ) 

I F( IL.EQ.O) GO TO 50 
WRITE(6.45) 

45 FORMAT(* ',4X, 'TEMPERATURE LIMITER HAS BEEN USED') 

IL=0 

50 IF( IF.Eq.O)GO TO 50 
WRI TE(6,55) 

55 FORMAT(* ' ,4X, 'CURRENT LIMITER IS IN USE') 

60 WRI TE(6,65) 

65 FORMAT( ’ ' ) 

WRITE(6,66) 

66 FORMAT( ' ' ,4X, 'TEMPERATURE PROFILE') 

WRI TE(6, 70) 

70 F 0 RMAT( ' ',9X,'0',11X,'r,11X,'2',11X,'3',11X,'4’,11X,'5',nx, 
C’6',11X,'7',11X,'8',nX, '9') 

WRI TE(6,65) 

00 100 1 = 1,91, 10 
IA=I-1 

WRI TE( 6, 80) I A, TA( I ), TA( 1 + 1 ), TA( M2) , TA( 1+3 ), TA( 1+4) , tA( 1+5) , 
CTA( 1+6), TA( 1+7), IA( l+8),TA( 1+9) 

80 FORMAT( ' ' , 1 2, 2X, F 10. 4 , 2X , F 10 . 4 , 2X , F 10 . 4 , 2X, F 10 . 4 , 2X, F 10. 4 , 2X, 
CF10.4,?X, F10.4.2X. F10.4,2X, F10.4,2X, F10.4) 

100 CONTINUE 
DTEMP^O.O 
WRITE(6,65) 

WRI TE(6, 1 10) 

110 FORMAI(' ',4X, 'CENTER TEMPERATURE AFTER EACH TIME 'NTERVAL') 

WRI TE(6, 120) 

120 FORMAT!' ' ,9X, ' I ' . 1 IX, '2' . 1 IX, ’ 3 ' . 1 IX. ' 4' , 1 IX. '5' , 1 IX, 
C'6'.11X.'7',11X,'8'.11X,'9'.11X,'10') 

WRITE(6,65) 

DO 140 1=1,91,10 
IA=I-1 

WRI TE(6, 130) lA, TC( I ) . TC( 1+1 ; . TC( I +2 ) , rC( 1+3 ) . TC( 1+4 ) , TC( 1+5 ) , 
CTC( l+6),TC( 1 + 7), TC( 1+8), TC( 1+9) 

1 30 FORMAT! ' ' . I2.2X. F10.4,2X. F10. 4.2X. n0.4.?X. F10.4.2X, F10.4.2X 
CF 10.4,2X. F10.4.2X, F10.4,2X. F10.4.2X, Flu. 4 ) 

140 CONTINUE 
WRI IF! 6, 65) 

WRI 1E!6, 150) 

150 FORMAI!' ',4X,' TOTAL CURRENT AFTER EACH TIME INTERVAL') 

WRI TE!6, 160) 

160 FORMAI!' ’ ,9X, ' 1 ' , 1 IX. '2' . 1 IX. ' 3' . 1 IX, '4' , 1 IX. '5' , 1 IX 
C'6',I1X,'7'.11X,'8'.11X.'9'.11X,'10') 

WRI TF!6.65) 

DO 180 1=1.91,10 
IA=I-1 

WRI TE!6. 170) lA. TT! I ) . 1 1! 1*1 ). TT( 1*2) . rr( 1 + 5) . TT! 1+4), TT! 1+5) . 
CTT| I+6),TT! 1+7), T’^ I 1+8). TT! 1+0) 
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170 FORMAT! ' \ I 2, 2X, MO. 4 , 2X . F 10 , 4. 2X, F 10 . 4 . 2X, F 1 0 , 4, 2X, F 1 0 . 4 , 2H, 
Cno.4,?X, no.4,2X, F10.4.2X, F10.4,2X, F10.4) 

180 CONTINUE 
WRITE(6.65) 

WRITE(6. 190) 

190 FORMAT! ' ',4X, ’VOLTAGE DURING EACH TIME INTERVAL’) 

WRI TE! 6,200) 

200 FORMAT!’ ’ . 9X/ 1 M IX, ’ 2 M IX/ 3 M 1X. ’ 4’ , 1 1 X, ’ 5 M 1X, 
C*6’,1lx//',11X/8’,11X/9’,nx/lO') 

WRITE!6,65) 

DO 220 (=1 ,91 , 10 

IA=I-1 

WRI TE! 6,210) I A, TV! I ), TV! 1 + 1 ), TV! 1+2), TV! 1+3 ), I V! 1+4 ). TV! I -►5), 
CTV( I +6 ) . TV! I +7 ) . TV( I +8 ) , TV( I +9 ) 

210 FORMAT! * ’ , I 2 , 2X , M 0 . 4 , 2X , f 1 0 . 4 . 2X, F 1 0 . 4 , 2X, F 1 0 . 4 , 2X, F 1 0 . 4 , 2X, 

CF10.4,2X, F10.4,2X, H0.4,2X, F10.4,2X, M0.4) 

220 CONT I NUE 

, 100,1 
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DO 230 1=1 , 
TC! I ) = 0 
TT( I )=0 
TV! > )=0 
CONTINUE 
RETURN 
FNO 

SUBROUTINE 


CURDEN 

Q ## ################### ### ^ 

###################### ;♦########## 
CALCULATION OF THE CURREN TDENS I TY 
0##############################^ #######^####### 

Q 4 ####### 4 ################ >>######### 


INTEGER I . I F, k. I ,M. I A, 1 1 . IM 
REAL TJ,WK,CD,DR, V. I TOT, TB. TA,VMAX,TC. TT 
C , 1 MAX ,DT,X,Al,Bt,Cl,Dl,nFTA. GAMMA , D E 
C,HL, IRR, TINF,WSI , PI ,C,A.UEL. LF, TV. V.,A 
COMMON I , I F, K, L,M, I A, 0, E. I L , A. IM 
C, TJ ,Wk,CO,OR. V. I TOT , TB( 101 ) , TA( 101 ) , VMAX, TC( 100) , TT( 100) 

C. IMAX, or ,X,A) ! in()),BI( lOUl.CI! 100), 01! 100), BETA! 100), GAMMA! 100) 
C,Hl , IRR. TIN! , WSI , PI ,C. DEL . I F. TV! 100) , VAR 
0- VAR + Lf *^EXP( - I 

COO .933E-2+8. 127F-6^TJ + 1 . 4*5 3E-7^EXP( 4 . 585E-2*TJ+0 ) 
C+2.297E-t?*tXP( 1 .20«E-2*rj^»U)*^V 
X- . 1 382*^V+.023 
CD^AMIN1!CD,X) 

RE TURN 
END 
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//RAHRTt JOB (0913-2-992-OU- . 33 : 00 . 900 ) , ‘ ^OB HARTMAN* 

// exEC FTG1CLG 
//C.SYSIN DO ^ 

######*#****#*#*##* #*-4*#*##** 

ENCAP ** 

C** THIS PROGRAM CALCULATES THE TEMPERATURE PROFILE OF A BACKWARD- ** 
C** BIASED SILICON SOLAR CELL. THAT IS ENCAPSULATED. ** 

C** A GAUSSIAN DISTURBANCE WITH DIAMETER DEL AND AMPLITUDE LF IS A ** 
C** PARAMETER THAT CAN FIT THE MODEL TO A MEASURED l-V CURVE. 

C** 

Q ♦***###*##**#*##*#*##■•»** 4 4-IMMf #* -IHHHf ** *4** ####*###»## <4 *### <4*# *# ##* * 

c** MULTIPLE RUN CONTROL ** 

^^*4#**4*4**4###*#4*#*4#*444#**44 444*4*4*4*44#*4#**444**4**#*4<»Mf4*444#4444 
Q**#4#4##»*444**#*4*#*4*4 44*#44*4444444444*#44#44*4 4#*#4 44*#4*#44##4*#4* 

INTEGER I . IF.K.L.M. lA. IL. IM 

REAL TJ.WKE.WKS.CO.DR, V. I TOT.VMAX, TC, TT . TE , AS. AE 
C.IBS, IBE. TAS. TAE.AIS.BIS.CIS.DIS.BETAS.OAMS 
C. IMAK.DT.X.AIE.BIE.CIE.DIE.BETAE.GAME.D. E 
C.HECE.HEEA, I RR . TA I R. WS. P I .C.A.DEL, LF .TV.WE, KE.VAR 
COMMON I, IF.K.L.H. lA. II , IM 

C.TJ,WKE.WkS.CO.OR.V, 1 TOT . VMAX. TC( TOO ) . TT( TOO ) . TE( 100 ) 

C.TBS( 101 ).TBE( 101 ) .TAS( 10 1 ) . TAE( 101 ) . AS . AE 

C.AIS( 100) .B!S( 100) ,CIS( 100) .DIS( 100 ) .BETAS( 100) ,GAMS( 100) 

C. IMAX.0T.X.AIE( 100 ) . B I E( 100 ) . C I E( 100 ) , D I E( 1 00 ) 

C,BETAE( 100),CAME( lUOl.D.E 

C.HECE, HEEA. IRR,TA1R,WS. P I , C. A. DEL. LF , TV( 100 ) . WE. KE, VAR 
C WE=.3B7 

C KE=7,24E-3 

WE=.4191 
KE=2.039E-3 
VAR=1 .C 
LF^n.6 
WE=WE/2 
CALL MN 
WE=WE^2 
VAR-1 .0 
LF=0.2 
CALL MN 
LF=0.6 
VAR=1 . 1 
CALL MN 


END 

SUBROUTINE MN 

(^##*#44#444#44# 4*44 #4#»4 444## 444# 4*4444##444#4#44*4**44*4*4#4*44#4*#44## 
Q4#44444*#444##44t44444#4 44444444 4444K 444444444444444444444 44444444444444 

MAIN 

Q4 44 444444 44444444444444 44444444 4 44444 4444444 4 44 4 4444444444444 4 444444 444 
Q4444 ♦444 44444 4444444444444 4444 4444444444444 4 4444 444 4444444 4444444444444 

INI EGER I . I F, K. L.M, I A, I [ , IM 

REAL TJ.WKE.WKS.CD.DR, V. I TOT . VMAX . TC , T'^.IE.AS.AF 
C. TBS. 1 BE. TAS. TAE . A I S . B I S . C 1 S , D I S , BE TAS . GAMS 
C. I.T\X. DT.X.AI E.BI E.CI E, Dl F. BETAE.CAME. D. E 
C.HECE.MEEA. I RR . T A I R . WS . P I . C , A . DEL . LF . TV. WE . KE . VAR 
COMMON I.IF.K.L.M.IA.IL.IM 

C, TJ.WkE.WKS.CD.OR.V, ITOT.VMAa. TC( U)0),TT( 100),TE( 100) 

C, TBS( 101 ) . TBE( 101 ) . TAS( 101 ) . TAf ( 101 ) , AS. AE 

C.AlSi 100).BIS( UU1),C1S( 100).DIS( l()0).Bf TAS( 100),GAMS( 100) 

C, I MAX, or ,X, Al E( 100 ) , B I F( 100 ) . C 1 E( 100 ) . D I E( 100 ) 

C.OETAEi 100 ), GAME ( 100),D,E 

C. HECE. HEEA, I RR. TA I R . WS. P I ,C. A.OEL. LF, TV( 100 ), WE. KE.VAR 
IM: 0 


DEL=23 

HECE=KE/( .3*WE) 
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o o r> o o n o n n n o n o o 


m rA^,o03 
I . \ 

TAIR'?*?.!) 

WS= . 04 

PI -3. l41^9Pr>54 
0 

ILnO 

DR=.0'?9602649 
C^lRR/HtLA+TAlR 
DO ^ ,76, 1 

TAS( I )-C 
TRS( 1 )rC 
^ CONI INUE 

00 in 1:^1 , 101 , 1 

lAf ( I )^C 
TBE( I )=C 
U) CONI INUE 
V 0.0 
VMAX~5»0.0 

CALCULAIION Of lOTAL CURRENl 


I A- 1 
I -0 

I TOT^O 
rj=TAS{ lA) 

CALCULATION OF CURRENTDENS I TY FOR 1=0 

CALL CURDFN 

MOT = (CD^PI^UK*^2)/4 

DO loo I ^ 1 . /'j, 1 

IA= 1+1 

TJ- TAS( 1 A) 

CALCULATION 0» CURKEN 1 DL NS I T Y FOR Kl^76 
CAl I CURDI N 

I TOT 1 10T+CD^I*P1*2^0R*^**P 
160 CONTINUE 

CM L Oil! PUT 
V' 2.00 
IMAX-- 1 . 

AS-2U 

AL-2 

1F=1 


RUNNINC CONTROL 

4f ## ff ## ff ff ## If ff 4f # ff ff ff ## 4f ff # ff 4f ###■** ff 4f ## <ff #### 4f # ff # ff 4f «f ff ff # fMf # 4f ## ff ff # ff «f <«f ### ff ff # ff ff If ff # # 


DO lao K- 1 ,40, 1 
I MAX- IMAX+.OS 
CAl L RUN 
CAl L RUN 
CAl I fINL 
CAL L OUl PU1 
lao CONT INUE 
R( lURN 
FNO 

SUBROUT INf RUN 

(^##ffff«ffff**f«*##H*fff*fH«Hf H«f# tfffH** IH»mf#ffH»f*fffff#<f«ftf*f4< tt<fffH#ff*f •♦♦Mf«f*fff*f#«fff«fM'«fff«f<f<fH# 

(^'‘«f##ff-if-H4fffffff#ff####-Hff4f«flfff#-^ff-Hff‘f*f4fMff#ffff4f4f<tf4f#-4ff««-HH4fff*f#ffffff«fff#fffffflfff#ffff#4fff#ff4fff4f 

raw ILMPLRAIURL RUN AND SURSUiNG 

^"#fflf4f4fH#ff-tfffffffffffff#«f«f^ffffff4f#^-if4f#«ftf#ff#####ff«f4fff-if-tfff#4fff#^-tfff##4f4f4f4f*f'«fff###H4f4f-tfff### 
(^^#ff##ffff#‘»f ff ff ff ff ## # #ff ff ff #f*<f ##ffffff *» # ###ff-ffff<f #ff#ffffff#ff ff4fff4f ff###1f-tf###ff ff ######### tf #ff 
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INTEGER I , I r . K. L .M. lA, I L. IM 

Rt aL TJ , WKE , WKS.CD. dr. V, I TOT , VHAX. TC . T T , TE . AS . AE 
C TBS. TDE. TAS. TAC . Al S. B I S, C I S. Dl S. BE T AS. CAMS 
c! IMAX.DT.X. AiE.BIE.CIE.OIE.BETAE.GAME.D, E 
C.HECE.MEEA, I RR. TA I R . WS, P I .C.A.DEL.LF, TV, WE, KE, VAR 
COMMON I, IE.k.L,M. lA. IL, IM 

C. TJ.WKE,WKS,CO,DR,V, I TOT.VMAX. TC( lUO) , TT( inO),TE( TOO) 
C,TBS( 101 ),TBE( inl ).TAS( 101 ).TAE( 101 j.AS.AE 
C Al S( 100 ) . Bl S( 100 ) .C I S( 100 ) , Dl S( lOU ) . BETAS( TOO) ,GAMS( 100) 
C. IMAX,DT,X,AIE( 1 00 ) . B I E ( 1 OO ) , C I E ( 1 00 ) , D I l ( 100) 

C.BETAE( 1 00). GAME ( 100).0,E 

C,HECE.HEEA, I RR. TA I R. WS, P I , C, A, DEL, LF , TV( 1 00 ) . WE , KE , VAR 
DO 161 L^l ,20. 1 
0T=.01 
CALL NEWIEM 

r oALL OUTPUT 

161 CONTINUE 

C CALL OUTPUT 

00 162 L=21 .40, 1 
0T= .005 
CALL NEWTEM 

C CALL OUTPUT 

162 CONTINUE 

C CALL OUTPUT 

DO 163 L = i:1.60, 1 

DT=.004 

CALL NEWTEM 

163 CONTINUE 

C CALL OUTPUT 

00 164 L=61 ,80. 1 

0T=.OO3 

CALL NEWTEM 

164 CONTINUE 

C CALL OUTPUT 

00 16^ L^81 . 100, 1 

DT=.0025 

CALL NfWTEM 

165 CONTINUE 

C CALI OUTPUT 

00 166 L=1 .20, 1 

DT=.002 

CALL NfWTEM 

166 CONTINUE 

C CAI L OUTPUT 

DO 167 L=21 .40, 1 
OTr.noi') 

CALL NfWTEM 

167 CONTINUE 

C CALL OUTPUT 

'30 168 L^41 .60. 1 

0T=.OOl 

CALL NEWTEM 

168 CONTINUE 

C CALL OUTPUT 

00 16Q L=61 ,80.1 

01^.001 

CALL NIWTLH 

169 CONI INUE 

C CALI OUTPUT 

DO 1 70 L:^81 . 100, 1 
DT- . 0008 
CALL NTWTEH 

170 CONTINUE 

C CALL OUTPUT 

DO 175 1=1 . 100, 1 
T1( I ) = 0 
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TC( ! )-0 
TV( I )^0 
175 CONTINUE 
C CALL OUTPUT 

DO 181 L=1 .20. 1 
DT=.00U6 
CALL NEWTEM 

181 CONTINUE 

C CALL OUTPUT 

DO 182 L=21.40.1 

0T=.0n04 

CALI NEWTEM 

182 CONTINUE 

C CALL OUTPUT 

DO 183 L=41 .60. 1 

DT=.0n02 

CALL NFWTEM 

183 CONTINUE 

C CALL OUTPUT 

DO 184 L=61.80.1 

dt=.ouoi 

CALL NEWTEM 

184 CONTINUE 

C CALI OUTPUT 

DO 185 L=81. 100.1 
DT= .onoi 
CALL NEWTEM 

185 CONTINUE 

C CALL OUTPUT 


RETURN 

END 

SUBROUTINE FINE 



INTEGER I . IF.K.L.M. lA, IL. IM 

REAl TJ.WkE.WKS.CD. DR. V. I TOT , VMAX. TC , T T , T E . AS , AE 
C, TBS. TBE. TAS. TAF , Al S. BIS. CIS. 01 S. BETAS. GAMS 
C. IMAX.DT.X, AIE.BIE.CIE.DIE.BETAE.GAME.D. E 
C, HECr . HEEA, IRR, TAI R,WS. PI . C. A, OFl . LF . TV. WE . KE. VAR 
COMMON I . I F. K. L.M, I A. 1 L. IH 

C TJ.WKr.WkS.CD.DR. V. I TOT .VMAX, TC( 100) , TT( TOO) , TE( 100) 

C. TBS( 101 ). TQE( im ) . TAS( U)1 ) . T AE ( 101 ) . AS , AE 
C A I S( 1 00 ) . B 1 S( UlO ) . C I S( 1 on ) , D I S( 1 00 ) , BE TAS( 100 ) , GAMS( 100 ) 
C. IMAX.DI .X.AIEI 100).BIE( 100).CIE( 10U).DIE( 100) 

C.BflAEI UlO). GAME ( 10O).D,E 

C. HFCE . HEEA. IRR. lAIR.WS. PI . C, A, DEL. LF, TV( 100 ) .WE. kE.VAR 
DO 175 1=1. 100.1 
TT{ I )=0 
TC( I ) = l) 

TV( I )‘0 
175 CONTINUE 

DO 181 L=1 .20, 1 

nr^ .00008 

CALL NEWiEM 

181 CONTINUr 

C CALL OUTPUT 

1)0 182 L^21 . 40. 1 
DT=. 00006 
CALL NEWTEM 

182 CONTINUE 

C CALL oulpur 

DO 18 i I -41 .60. 1 
0^ 00005 
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oonoooon C300000 00000 


CALL NLWTEM 

183 CONTINUC 

C CALL OUTPUT 

DO 184 L=61 .80, 1 
0T=. 00003 
CALL NLWTEM 

184 CONTINUE 

C CALL OUTPUT 

00 185 L^81 , 100. 1 
0T =. 00001 
call NfWTEM 

185 CONTINUE 

C CALL OUTPUT 

RETURN 
END 

SUBROUTINE NEWTEM 


NEW TEMPERATURE CALCULATION ** 


INTEGER I . I r . K. I .M. lA. IL. IM 

REAL TJ.WKE.WKS.CO.OR.V. I TOT . VMAX. TC . TT , TE. AS.AE 
C. TBS. TBE. TAS. TAE . A I S , B I S , C 1 S. D I S . BE TAS . GAMS 
C. IMAX.OT.X.AIE.BIE.CIE.DIE.BETAE.GAME.D.E 
C.HfCE.MEEA. IRR. TAIR.WS. PI ,C. A. DEL . LF . TV. WE , KE . VAR 
COMMON l.lF.k.L.M.lA.IL.IM 

C. TJ.WKE.WKS.CD.OR.V. I TOT .VMAX. TC( 100 ) . TT( 100 ) . TE( 100 ) 

C. TBS( 101 ) . TBE( 101 ) . TAS( 10 1 ) . TAE( 1 01 ) . AS. AE 
C.AIS( 100).B1S{ 10O),CIS( 100).DIS( 100).BETAS( 100).CAMS( 100) 

C. 1MAX,DT.X.AIE( 100).BIE( 100),CIM 100).0IE( 100) 

C.BETAE( 100 ). GAME ( 100).D.E 

C.HECE.HECA. IRR, TAIR. WS. pi .C. a. del. LF, TV( 100 ) .WE. KE. VAR 
10 FORMAT!' *.'A \F10.4,' B \F10,4.’ C '.Flu. 4.' 0 *.F10,4) 

1=0 

IA=I>1 
TJ= TBS( lA) 

TK=TBC( lA) 

CALCUIATION OF CURRENTOENS I TY FOR 1=0 
CALL CURDEN 

calculation or thermal conduction 

WKS-WS/! ,58*^TJ/350) 

WKF=WF^KF 

calculation of THERMAL DIFFUSION 

AS=WKS/( 1 .647*WS) 

AE=WKE/( 1 ,28^WE) 

CALCULATION OF FIRST MAIRIX ELEMENTS ON THE CELL 
AIS( 1A)=0 

BIS( lA) = -?/DK»»^2-1/!AS*»Dl ) 

CIS! IA)-,VOR’*^:? 

DIS( IA)-(2/DR*^^P-l/! AS^Ol ) ) ♦► TJ T BS! IA-»*1 )/DR^^2 
C-iCD^V-HLCC^^Tj^HECE’^TKl/WKS 

C CALCULATION OF THE FIRST MATRIX ELEMENTS ON THE ENCAPSULANT 

C 

Al E( lA) 0 

BIE( IA) = -?/0R^*^2-1/(AE*DT) 

CIE! IA) = P/DK**^P 


:35 









c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


DIE( IA) = (3/DR^*2-1/( AE»»DT) )^TK-2<^TBE{ IA+1 )/DR**2 
C-(HECE*TJ-HECE*^TK-HEEA*rK-^HEEA^TAIR+lRR)/WKE 
WRITE(6, U))A1S( IA).BIS( IA).CIS( IA).DIS( lA) 

WRITE(6, 10)AIE( IA),BIE( IA),CIE( IA),0IE( lA) 

00 100 1=1 . 7U, 1 
IA=I+1 

TJ = TBS( I A) 

TK=TBE( IA) 

CALCULATION OF CURRENTDENS I TY FOR 1=1 TO 1=74 
CALL CUROEN 

CALCULATION OF THERMAL CONDUCTION 
WKS=WS/( .58+TJ/330) 

CALCULATION OF THERMAL DIFFUSION 
AS=WKS/( 1 .647^WS) 

CALCULATION OF THE MATRIX ELEMENTS ON THE CELL 

AIS( IA)=( 1-1 )/( (2*1-1 )*0R**2) 

BIS( IA) = -1/0R**2-l/( AS^^DT) 

C1S( IA;= l/( (2^1-1 )*DR*^2) 

0IS( I A) = -AIS( IA)^T0S( IA-1 ) + ( 1/DR^^^2-1/(AS^^0T) )^^TJ 
C-CIS( IA)*TBS( IA*^1 )-(CD*V-HECE*TJ+HECE*U.)/WkS 
C+( TBS( I A+1 )-TBS( IA-1 ) )*^2/( ( 205 + TJ )<^4*0R**2) 

CALCULATION OF THE MATRIX ELEMENTS ON THE ENCAPSULANT 

AIE(IA)=(I-1)/((2*I-1 )*0R**2) 

BIE( IA) = -1/0R^*2-1/( AE^^OT) 

CIE( IA)=l/( (2*1-1 )*nR**2) 

0IE( IA) = -A1E( IA)*TBE( I A- 1 ) + ( 1 /0R**2- 1 / ( AE*OT ) ) *TK 
C-C1E( I A)*TBE( IA+1 )-( +HECE*Tj-HrCE*TK-MtEA*1k+HEEA*TAIR+lRR)/WKE 
WRlTE{6. 10)AIS( IA).BIS( IA).CIS( IA),DlS( IA) 

WRI TE(6, 10)AIE( IA).BI£( IA).CIE( IA).0IE( IA) 

100 CONTINUE 

1 = 7 ^ 

IA=I+1 
TJ=TBS( IA) 

Tk=TBE( I A) 

CALCULATION OF CURRENTDENS I TY FOR 1=75 
CALL CURDEN 

CALCULATION Of THERMAL CONDUCTION 
WKS=WS/( .38+TJ/3SH) 

CALCULATION OF THERMAL DIFFUSION 
ASrWkS/( 1 .f>47*WS) 

CALCULATION OF MATRIX ELEMENT FOR 1=73 ON THE CELL 

A1S( I A)'--1/(2*dR**2) 

BIS( IA)r..T/(?*DR**2).i/(AS*DT) 

C I S( I A) -n 

0IS( 1A) = -AIS( 1A)*TBS( IA-n-(BIS( ! A ) +2/ ( AS*OT ) ) *T J 
C-(CD*V-HLCL* IJ + HECE'^T k )/WKS 

C4-73. 3*MErE*WS*( 1 .3*TBS( I A)-, 5'*TBS( I A-1 )-TRE( I A+1 ) ) / ( 74 . C*WkS*DR ) 
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c 

c 


c 

c 

c 


c 

c 

c 


c CALCULATION OF MAff^lX ELEHtNIS FOR 1 = 75 ON THE .NCAPSULANT 

AIE(IA)=(I- 1 )/(( 2 »|-i )*DR** 2 ) 

Bi£( IA) = -1/0R**2-1 /( aE^PT) 

CIE( IA)=|/( (2*1-1 )*0R**2) 

= IA)*T8E( IA-1 ) + ( 1/DR**2-1/( AE»DT1 l*Tk 

RRirfi«,io)Ai-iiA}.BiE^A!:ciE(i;!:6l!liA| 

IA=I+1 

CALCULATION MATRIX ELEMENT FOR 1=76 ON THE ENCAPSULANT 
AIE( IA) = ( l-1)/((2*l-1 )*nRi»#2) 

B1 E( IA) = -1/DR**?-1/( AE^^DI ) 

CIF{ IA)=I/(( 2 *|-i 

Dl E( I A)=-AI E( IA)^TBE( IA- 1 )+( 1 /DR^* 2 - 1 /( AE^DTl \a\ 

C-CIE( IA)*TBE( IA+1)-(-HEEA*TBE( lAl+HEFAiTAIRl RRw5ir ’ 

I A= I + 1 

CALCULATION MATRIX ELEMENTS ON THE ENCAPSULANT 

AIE{IA) = (|-1)/((2#|.-| 

BIF( IA) = -l/OR^*2-l/( AE^OD 
CIE( IA)=|/((?*|-i )»DR*^2 ) 

c ^WR >A)fHEEA^TAIR+ RR)/WKE^ 

150 CONtUuE ’ ‘ IA).DIE( lA) 

j CALCULATION OF TEMPERATURE ON THE CELL AT THE END OF PERIOD 
I A- I 

BETAS( IA)=B1S( lA) 

CAMS( IA)=OIS( IA)/BIS( lA) 

DO 380 IA= 2 . 76.1 

IA) = BIS( IA)-(AIS( IA)*CIS( IA-1 ))/BE TAS( I A- 1 1 

TAS( 7Tj ) = GAMS( 76 ) 

DO 43U M-1 , 1 

I A = 76-M 

U30 CONIINUE^^”^* IA)-(CIS( IA)HTAS( IA+1 ) )/BETAS( lA) 

CALCULATION OF TEMPERATURE IN THE ENCAPSULANT AT THE END 01 PERIOD 
lA^l 

BETAE( IA)=0IF( IA) 

CAME( lA) = On:( IA)/BIE( IA) 

DO I A- 2. ino. I 

IA)-BIF( IA)-(AIE( IA)^CIF( IA-1 ll/BETAFl IA-11 

435 CONTIn0e"‘°'^‘ IA)-AIL( IA)*CAME( IA-1))/BLrAL(IA) 

iA=inn-M 

436 CONtInuE^^'^^'^ ' 'A)-(CIE( IA)*TAE( IA-M ) )/BETAE( IA) 

TAE( 101 )=TAF1 I"-’; 

temperature limitlr 
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DO ^^0 1^1,76,1 

IF(TAS( I I.LT.IOOO) GO TO 440 

\L--} 

TAS{ \ )=1000 
440 C0NHNUE 


CALCULATION OF TOIAL CURRENT 




#«###«####################«#«#»## 


1=0 

IT0T=0 
TJ=TAS( 1 ) 

CALCULATION OF CURREN FDENS I TY FOR 1=0 

CALL CURDEN 
lT0T=(CD*PI*»DR^*2)/4 
DO 430 1=1,75,1 
IA= 1+1 
TJ=rAS( lA) 

CALCULATION OF CURRENTDENS I TY FOR Kl<76 
CALL CURDEN 

I ror=i Tor+CD*i^Pi*2^0R**2 
430 CONTINUE 

TC(L)=TAS( 1 ) 

TE{L)=TAE(1) 

Tr(L)=lTOT 
TV( L)=V 

CURRENT LIMITER 

##♦###########*#### 

c 

D= 1 . 0000000 1 ^ ( MAX 
99999999# I max 
IF ( I TOT.GT. D) GO TO 490 
IF { I ror.LT.E.AND. IF.EQ. 1 ) GO TO 495 
GO TO 300 

490 V=(V*IMAX)/I TOT 
I F = 1 

GO TO 300 

493 V={V^IMAX)/I TOT 

I F( V, LT. VMAX) GO TO 500 
V^VMAX 
I F -'O 

300 CONMNUE 
IM= IM+1 

DO 320 1=1 , 76, 1 
TBS( I ) = TAS{ I I 
520 CONTINUE 

DO 340 1 = 1 , 101 . 1 
TBE( I )-TAE( I ) 

340 CONI INUE 
C CAU OUlPUT 
RETURN 
END 

SUBROUTINE OUTPUT 

(^*i‘##-*########-tf######*< #*♦######■•♦ ##############################-»*#### »■###### 

(^####### ######################### #4# ###########■##### 4^# 4####+*# 

ouirui IN rAiu.L iorm 

Q##4#######<############4##4##4###4#####4######################444####4## 

( j ######'»#####- H '#»-#« f #### i » 4 ## #♦>##’»♦############ ######################## 

INTEGER I . I F. K, I..M. 1 A. I I . IM 

REAL TJ.WKE. WKS.CD, DR, V. I TOT. VMAX, T C , I T . I E . AS , AL 
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C, TBS, lOE.TAS. TAt,AIS,BIS.CIS,OIS,BEIAS.CAMS 
C. IMAX.nT.X.AIC,Bir,CI£.OIE,BETAE,GAME,D,E 
C.HECE.HEEA. I RR , 1 Al R, WS, PI . C, A. DEL. LE. T V, WE . HE, VAR 
COMMON I . lE.h.l .M. lA. IL. IM 

C, TJ.WkE.WKS.CO.OR.V. I TOT, VMAX. TC( 100).TT( 100). TE( 100) 

C, TBS( 101 ),TBE( 101 ), TAS( 101 ).TAE( 101 ).AS,AE 

C, AIS( 100),BIS( 100),CIS( 100),DIS( 100).BETAS( 100),CAMS( 100) 

C, IMAX.OT,X,AIC( 100).BIE( 100),CIE( 100),0IE( 100) 

C,BETAE( 100), GAME ( 100).0,E 

C,HECE,HEEA, IRR,TAIR,WS, PI , C, A, DEL, LF, TV( 100), WE. KE,VAR 
C 

C PRINT STATUS 

C 

IF( IM.EQ.O) GO TO 1 
GO TO 9 

1 WRITE(6,2)IM,VAR 

2 FORMAT( ' 1 ' ,'4X, ' I IME = ’.I12.2X,' INCREMENTS, THE VARIATION', 

C' FROM THE STANDARD CURRENT IS ’,F12.4) 

WRITE(6.3jDEL.LF 

3 FORMAT!' '.ilX.'THE GAUSSIAN IS ',F12.i|. 

C WIDE AND HAS AN AMPLITUDE OF ',F12.4) 

CO TO 15 
9 WRI TE(6, 10) IM 

10 FORMAT! 'I'.UX,' TIME = ',I12,2X,' INCREMENTS') 

15 WRITE!6,20)V 

20 lOfMAT!' ' ,4X, 'VOLTAGE = ' , F 12. 4 , ' V. ' ) 

WRITE!6,30)IT0T 

30 FORMAT!' ' ,4X. 'CURRENT = ' , F12. 4, ' A. ' ) 

IF! IL.EQ.O) CO TO 50 
WRITE!6,45) 

45 FORMAT!' ' ,4X. 'TEMPERATURE LIMITER HAS BEEN USED') 

IL=0 

50 I F! I F. EQ.O)GO TO 60 
WRITE!6,55) 

55 FORMAT!* ' .4X. 'CURRENT LIMITER IS IN USE') 

60 WRITE! 6, 65) 

65 FOR.MAT! ' ' ) 

r 

C PRINT TEMPERATURE PROFILE OF THE CELL 
C 

WRI TE! 6. 66) 

66 FORMAT!' ' .4X. ' TEMPERATURE PROF ILF OF THE CEIL') 

WRITE!6. 70) 

70 FORMAT! ' '.9X.'0',11X.'r.11X,'2'.11X,'3'.11X.'4'.11X.'5'.11X, 

c'6' . 1 IX, '7' , 1 1X, '8' , nx. '9' ) 

00 100 1=1.91,10 
IA-l-1 

WRI It !6,80) I A. TAS( I ) . TAS! 1*1 ) . TAS! I +2 ) . TAS! I +3 ) , IAS! I +4 ) . TAS! I +5 ) , 
CTAS! 1+6), IAS! !♦/). IAS! 1+8), IAS! 1+9) 

80 FORMAT! ' ' , 1 2 . 2X. M 0 . 4 . 2X. F 1 0 . 4 , 2X , F 10 . 4 , 2X, F 10. 4 . 2X, F 10. 4 , 2X, 

CF10.4.2X, F 10.4, 2X. F 10.4. 2X. F 10.4, 2X, FI 0.4) 

100 CONTINUE 
WRI TE!6.65) 

C 

C PRINT TEMPERATURE PROFILE OF THE ENCAPSUI.ANT 

C 

WRI TE(6, 101 ) 

101 FORMAT!' ' ,4X. 'TEMPERATURE PROFILE OF THl ENCAPSULAN T ' ) 

WRI IE!6, 102) 

102 FORMAI!' '.OX.'o'.lIX ' 1 ' . 1 IX, '2' , 1 IX. ' 3' . 1 IX, '4' , 1 IX, '5' . 1 IX. 

C'(j' . 1 IX, ' 7' . 1 1X. '8' , 1 IX, '9' ) 

DO 10.1 1 = 1.91,10 
IA=I-1 

WRI FF! 6. 104 ) I A. lAF! I ) . 1A1 ! I + 1 ) . IaF! I♦■'’ ) . lAF! 1 + 3 ) . lAF! I +4 ) , 

CTAEI I+5),TAE! 1+6). lAE! 1+7,, lAEI 1+8). lAE! 1+9) 

1( 4 FORMAT! ' ' . I2.2X, FI0.4..?X, n0.4.2X, F10.4.2X. F10.4.2X, FU). 4.2X, 
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Cnn.i4.?K. F10.4.2X. F10.4.2X, nO.U,2X, M0.4) 

103 CONI iNUe 
WRITE(6»65) 

H. STORY OF CENTER TEMPERATURE OF THE CELL 
WRITE(6. 110) 

no FORMAT^* \4X/CENTFR TEMPERATURE OF CELL * 

C. ’AFTER EACH TIME INTERVAL*) 

WR1TE(6. 120) 

1 20 FORMAT! ' ’.9X,'l'ni>^.*2',nX,*3’,11X/4*.nx.*5*.11X. 

C'6* nix, *7' mx. *8 \ 11X, '9M1X, ' 10' ) 

DO 140 1=1,91.10 
I A=l-1 

WRI TE(6. 130) I A. TC( I ) . rC( I +1 ) . TC( n2) nC( 1+3 ) . TC( 1+4) . TCU+5 ) , 
CTC( 1+6). TC( l+7).TC( l+8).TC{ 1+9) 

130 FORMAT! ' ’ . I 2, 2X . F 1 0 , 4 , 2X . FI 0 . 4 , 2X , F 1 0 . 4 . 2X, F 1 0 . 4 , 2X, F 10 . 4 , 2X . 

CF10.4.2X. F10.4.2X, F10,4,2X. F10.4.2X. F10.4) 

140 CCNl INUE 
WRITE!6.65) 

HISTORY OF THE CENTER TEMPERATURE ON ENCAPSULANT 
WRI TE! 6, 141 ) 

141 FORMAT!* ’.4X, 'CENTER TEMPERATURE OF ENCAPSULANT ' 

C. 'AFTER EACH TIME INTERVAL*) 

WRI TE!6. 142) 

142 FORMAT ( ' ' . 9X . * 1 ' .1 1 X 2 ' ,1 1 ^ ' 3* . 1 1 X , ' 4 * ,1 1 X , * 5 * . 1 1 X. 

C'6’niXn7,11X.'8' .1ixn9* nix. *10* ) 

DO 144 1-1.91 . 10 
I A- 1-1 

WRI TE!6. 143) lA. TE( ! ) ,Tf ( I +1 ) .TE( 1+2) nE( 1+3 ) .TE( 1+4), TE( 1+5) . 
CTF( I +6) ,TE ! I + n . TE ( 1+8 ) , TE( I +9 ) 

143 FORMAT! ’ * . I 2 . 2X , M U . 4 . 2X . F 1 0 . 4 . 2X . F 1 0 . 4 , 2X. F ^ 0 . 4 . 2X. F 1 0 . 4 . 2X. 
Cr 10, 4.2X. F U). 4.2X, F 10. 4 . ;?X. F 10. 4 , 2X. F 10. 4 ) 

144 CONI INUE 
WRI TE! 6.65) 

HISTORY OF TOTAL CURRENT 
WR; TL!6. 150) 

150 FORMAT!* ’.4X.' TOTAL CURIUNT AITER EACH TIME INTERVAL*) 

WRI TE! 6. 160) 

160 TORMAl! * * . QX ,* .1 1 X . ’ 2 * . 1 1 X . ' 3 * .1 1 X 4 ' .1 1 X , ' 5 * . 11 X . 

C'6’niX.’7'.11X,'8'.11xn9'.nX.*10' ) 

WRI IE! 6.65) 

DU 180 1=1,91,10 
I A-- 1-1 

WR I I E! 6. WO) I A. TT ! 1 ) , T T ( I + 1 ) , IT ! 1 +2 ) . T 14 I + 3 ) , TT ! I +4 ) . TT! I +5 ) . 
C1T( t+6),TT! 1+7), TT! 1+8), T 1 ( 1+9) 

1 /n I ORMAT! ' ' , 12. 2H. I 1 O . 4 . 2X . f 1 0 . 4 . 2X , F 1 0 . 4 . 2X . f 1 0 . 4 . 2X , F 1 0 . 4 , 2X . 

CF 10.4.2X, F10.4.2X. FlU,4.2X, F10.4,2X, F10.4 ) 

180 CONTINUE 
C WRITE! 6.65) 

C HISTORY OF VOLTAGE 

C 

C WR I TF! 6. ion) 

190 I0RMAI(* ' .4X. 'VOLTAGE DURING EACH T I ML INTERVAL') 

G WR 1 Tl ( 6.200 ) 

200 fOKMAl( ' ’ .9X. ' 1' . in. ' 2 ' , 1 1\. n* . 1 1X.*4’ , 1 IX, '5' . nx. 

G ' 6 ' . 11 X , ‘ 7 ’ . 11 X , ’ 8 ’ . 1 1 X . ’ 9 * n 1 X , ' U) ’ ) 

C WRI TE 16.65) 

DO 220 1-1.91,10 
1 A- I -1 

WRI TE(6.210) 1 A. TV( I ) . rV( I + 1 ) . TV( 1+2) nV( 1+3 ) . TV! 1+4) . TV( 1+5 ) , 




C 



C crv( l+6).TV( l+7).TV( 1+8). rv{ 1+9) 

210 FORMAf( ' ' . I2.2X. ri0.4,?X, F10.4.2X. rU).4.2X, F10.4,2X. F10.4.2X, 
CF10.4.2X, F10.4.2X. F10.4,2X, F10.4.2X, F10.4) 

220 CONTINUE 

DO 230 1 = 1 , TOO. 1 
TC( I )=0 
TE( I ) = 0 
TT( I )=() 

TV( I )=0 
230 CONTINUE 
RETURN 
END 

SUBROUTINE CURDEN 

Q######*######### »#############»#####*######*#*######*######## )♦##*###### 

##*#*###### #######*+♦ ###*#**######*####»#*#*## 

CALCULATION OF THE CURREN TDCNS I TY 

########>> 'If 

lf#fl 4f#*###4f####-lf#<»f'lf#/Hf##'lf##'»f-lf*##<H If##*##* 

INTEGER I . I F. K, I ,M, lA, I L, »M 

REAL TJ.WKE.WKS.CD.DR.V. I TOT, VMAX, TC . TT , TE . AS, AE 
C, TBS. TBE, TAS, TAE . A I 3, B I S . C I S . D I S . BETAS. GAMS 
C. IMAX,DT,X,AlE.BIE.Ct E.DIE.BETAE.GAME.D, E 
C. HECE, HEEA, I RR . TA I R . WS, P I . C , A, DEL . LF , TV, WE , KE, VAR 
COMMON I , I F. K. 1. ,M. I A, I L. IM 

C. TJ ,WKE. WKS.CD, OR, V. I TOT . VMAX, TC( 100 ) , TT( TOO ) , TE( 100 ) 

C, TBS( 101 ) , TBE( 101 ) . TAS( 1 U 1 ) , T AE ( 1 01 ) . AS , AE 

C.AI S( 100) .BIS( 100) ,CIS( 100) ,0IS( 100 ) , BETAS( ’00) ,GAMS( 100) 

C. IMAX,DT,X,A)F( 100).BIE( 100),C1E( 100),DIE( lOO) 

C.BEFAEf 100).GAME( 100), O.E 

C. HECE, HEEA, I RR . TA I R , WS , P I . C , A , DEL , LF , T V( 1 00 ) . WE , KE . VAR 
D=VAR + Lr^EXP( - l^^2/DEL*^^2) 

C0= 1 .933E-2+8. 1 2 ?E-6*^T J + 1 . 453E- 7’*EXP( 4 . 585E-2+TJ ) 

C+2.297E-5^EXP{ 1 . 208E-2*TJ^O ) <^V 
X^. 1382^V+.023 
C0=AMIN1(C0,X) 

RETURN 

END 
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